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ABSTRACT 


The Lorentz force equation governing the interaction between radiation and 
charged particles in the presence of a static magnetic field can be reduced to the form 
of the simple pendulum equation, when applied to the Free Electron Laser (FEL) and to 
the amplification of whistler waves in the magnetosphere. 

The first topic in this thesis is the start-up of an FEL with a compact design having 
a small, spatial undulator period, many penods, and a small undulator field strength. 
Under these conditions, the number of photons spontaneously emitted into the 
resonator mode each pass is so small that a classical field cannot be established. 
Quantum fluctuations affect the wave-particle interaction and therefore diminish weak- 
field gain. An FEL start-up condition, several orders of magnitude more restrictive than 
determined by quantum recoil alone, is derived from the photon rate equation. 

The next topic is the evolution of a weak, classical radiation field dnven by a long 
electron pulse. Assuming low gain, it is possible to reduce the coupled optical wave 
and electron pendulum equation to a simple, first-order, partial differential equation 
including electron dynamics, slippage, desynchronism, and resonator loss. This 
analytical method is used to study coherence development, and the effects of varying 
the electron beam pulse shape. 

Next, the Boeing Average Power Laser Experiment (APLE) is investigated using 
two and three dimensional numerical simulations. Vanous methods of reducing the 
electron-radiation coupling in the oscillator are examined in order to minimize the 
electron beam energy spread introduced into the FEL amplifier. An alternative for 
achieving high power makes use of the oscillator alone, but with a tapered undulator to 
increase the electron-radiation coupling. 

The final topic is charged particles, following helical trajectones, traveling in the 
magnetosphere between mirror points along the earth's magnetic field lines. When 
electrons interact with a whistler wave near the geomagnetic equator, their motion can 
again be described by the simple pendulum equation. The inhomogeneity of the 
earth's magnetic field near the geomagnetic equator is analogous to the field in a 
tapered undulator. Numerical phase-space simulations are developed to investigate 


the interaction. 


iv 








Public reporting burden tor this collection of information 1s estimated to average | Nour per response, including the CMONTEREY: SIT sources, 
gathering and maintaining the data needed. and completing and reviewing the collection of information Send comme d ct of this 
eote ceo ohitlommashon, inidluding tuggenho for reduc Thi Cul sen lO WaEhUnglor Aes Qu ares Services, Directorate for int semston Operscrm ad Reports, 1215 jeherson 
Davis Highway. Suite 1204, Arlington: VA 22202-4302. and to the Office of Management and Budget. Paperwork Reduction Project (0704-0 188), Washington, OC 20503. 


1. AGENCY USE ONLY (Leave blank) | 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 
June 1993 Dissertation 
4. TITLE AND SUBTITLE : S. FUNDING NUMBERS — 


A Study of the Amplification of Laser and VLF 
Waves Using a Simple Pendulum Model 


6. AUTHOR(S) 


Keith Alan Sturgess 





7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION 


REPORY NUMBER 





Naval Postgraduate School 
Department of Physics 
Monterey, California 93943 


9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING / MONITORING 
AGENCY REPORT NUMBER 


11. SUPPLEMENTARY NOTES 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 1i2b. DISTRIBUTION CODE 


Approved for Public Release; distribution is unlimited. 





13. ABSTRACT (Maximum 200 words) mye Lorentz force equation governing the interactio 
between radiation and charged particles in the presence of a static 
magnetic field can be reduced to the form of the simple pendulum 
equation, when applied to the Free Electron Laser (FEL) and to the 
amplification of whistler waves in the magnetosphere. 

The first topic in this thesis is the start-up of an FEL with a 
compact design having a small, spatial undulator period, many periods, 
and a small undulator field strength. Under these conditions, the 
number of photons spontaneously emitted into the resonator mode each 
pass is so small that a classical field cannot be established. Quantum 
fluctuations affect the wave-particle interaction and therefore diminis 
weak field gain. An FEL start-up condition, several orders of magnitude 
more restrictive than determined by quantum recoil alone, is derived 
from the photon rate equation. 

The next topic is the evolution of a weak, classical radiation 
field driven by a long electron pulse. Assuming low gain, it is 


ieee hlectron Laser, FEL, undulator, start-up, oe GEECODE 
quantum me uatiens a meget modes, yi er or 


. SECURITY CLASSIFICATION 8. “SECURITY CLASSIFICATION : SECURITY CLASSIFICATION 20. LIMITATION OF ABSTRACT 
OF REPORT OF THIS PAGE OF ABSTRACT 


UNCLASSIFIED UNCLASSIFIED 
NSN 7540-01-280-5500 Standard Form 298 (Rev 2-89) 


aa se 2 = => 













UNCLASSIFIED UL 





REPORT DOCUMENTATION PAGE, SF 298 


igs 


Continued 


possible to reduce the coupled optical wave and electron 
pendulum equation to a simple, first-order, partial 
differential equation including electron dynamics, slippage, 
desynchronism, and resonator losses. This analytical 
method is used to study coherence development, and the 
effects of varying the electron beam pulse shape. 

Next, the Boeing Average Power Laser Experiment (APLE) 
is investigated using two and three dimensional numerical 
Simulations. Various methods of reducing the electron- 
radiation coupling in the oscillator are examined in order 
to minimize the electron beam energy spread introduced 
into the FEL amplifier. An alternative for achieving high 
power makes use of the oscillator alone, but with a tapered 
undulator to increase the electron-radiation coupling. 

The final topic is charges particles, following helical 
trajectories, traveling in the magnetosphere between mirror 
points along the earth's magnetic field lines. When electrons 
interact with a whistler wave near the geomagnetic equator, 
their motion can again be described by the simple pendulum 
equation. The inhomogeneity of the earth's magnetic field 
near the geomagnetic equator is analogous to the field ina 
tapered undulator. Numerical phase-space simulations are 
developed to investigate the interaction. 
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l. INTRODUCTION 


A free electron laser (FEL) oscillator uses a beam of relativistic electrons 
passing through a periodic undulator to produce and amplify light stored within a 
resonator. Since this light is a product of the electron - undulator interaction, its 
properties (i.e. wavelength and power) can be changed by altering the properties 
of the electron beam and the undulator. The ease of tunability of the FEL over a 
wide range of wavelengths is beginning to make the FEL a sought after piece of 
equipment by many disciplines in universities, research labs, and industry. The 
promise of high efficiency and high average power has made the FEL a 
research priority of the armed forces as they seek a speed-of-light weapon. 

The basic physics of an FEL is presented in Chapter Il as a starting point for 
the research chapters that follow. The FEL theory presented here uses 
Maxwells’ field equations to describe the evolution of the light in response to the 
bunching electrons, and self-consistently uses the pendulum equation to 
describe the evolution of the bunching electrons in response to the light [1]. 

Chapter Ill examines FEL physics during the first several passes of the 
electrons through a compact undulator that uses many short periods and a small 
undulator parameter. Spontaneous emission is introduced, and the effects of a 
small number of photons in the resonator causes statistical fluctuations of the 
photon number. These quantum fluctuations, if large enough, can significantly 
degrade the classical, small-signal gain of the FEL causing slow start-up or 
possibly the failure of the FEL to start. 

New and original research concerning the problem of quantum limitations on 
the start-up of FELs is presented in Chapter Ill. The previously derived equation 


for gain including photon fluctuations is incorporated into the photon rate 


equation to obtain a new start-up condition. The photon rate equation is then 
solved analytically using the new expression for gain. This and other analytical 
results are compared with those obtained by other authors using a completely 
quantum mechanical method. Also, the new expression for gain is incorporated 
into a numerical simulation to solve for optical mode evolution. 

Chapter IV follows the evolution of long optical pulses using a partial- 
differential equation which includes effects of desynchronism and resonator 
losses. Desynchronism describes the misalignment between the optical and 
electron pulses on Successive passes through the resonator. Analytical 
expressions are derived for the optical wave evolution assuming a long electron 
pulse, while the FEL operates in the small signal, low gain regime. The first- 
order, partial-differential equation is then used to examine the effects of the 
electron pulse shape on the final optical pulse. Also, the coherence 
development from a random incoherent optical pulse to a final coherent optical 
pulse is examined. 

A partial-differential equation describing optical pulse evolution is used for 
new and original research in Chapter IV. The single-pass, first-order partial- 
differential equation is extended from one to multiple passes including 
desynchronism and resonator losses. The partial-differential equation is solved 
analytically for the case of a continuous electron beam, a parabolic electron 
pulse shape, and a Gaussian electron pulse. A desynchronism curve showing 
the variation of optical power with desynchronism is then created from the 
analytical solution for a parabolic electron pulse shape and compared to 
desynchronism curves produced numerically, experimentally, and with analytical 
methods used by other authors. It is shown that exact synchronism gives zero 
optical power in steady-state. The case of a triangular electron pulse is solved 


analytically. A numerical FEL simulation is modified to solve the single-pass 


partial-differential equation and used to study coherence evolution. Finally, a 
numerical code is developed to compare the analytical solution and the 
numerical FEL simulation for the cases of a parabolic, triangular or continuous 
electron pulse. 

Only weak optical fields are studied through Chapter IV. Chapter V 
examines the Boeing/LANL high Average Power Laser Experiment (APLE). 
Since the proposed design of APLE has a single electron beam powering both 
an oscillator and amplifier FEL in series, it is important that the oscillator FEL not 
induce a large energy spread in the electron beam prior to the amplifier FEL. 
Chapter V examines methods of reducing the electron - light coupling in the 
oscillator. One method of reducing the coupling calls for inducing a ramp, or a 
slew, in energy along the electron pulse. Running the oscillator alone as a high 
average power FEL is also examined. 

The investigation into methods of reducing the electron beam energy spread 
induced by the FEL oscillator at the Boeing APLE project is new and original 
research. A numerical FEL simulation is extended to accurately calculate the 
optical field evolution in the presence of large desynchronism. Using this 
modification, a desynchronism curve is presented for the APLE oscillator. The 
desynchronism that best reduces the electron energy spread for the APLE 
oscillator is found. The simulation is modified further to include the effects of 
energy and timing jitter of the electron pulse from the accelerator. These effects 
are explored for large and small values of desynchronism. The simulation is 
modified again to include an energy slew across the electron micropulse in order 
to show that the energy spread induced in the electron beam is more than 
Boeing researchers had expected. The final APLE topic explores the tapered 
oscillator design as a high average power device. Conditions are obtained for 


achieving the optimum power from the FEL oscillator alone with no amplifier. 


Finally, Chapter VI examines whistler wave amplification employing the 
similar theoretical and simulation methods used to understand the FEL. The 
whistler and FEL both have electrons in helical trajectories coupling to and 
amplifying a transverse electromagnetic wave. The wave and electron 
equations of motion are developed for the whistler interaction including the 
inhomogeneity due to the variation in magnetic field strength along a field line of 
the earth's magnetic field. Simulations are then developed from these equations 
in order to describe stimulated VLF emission. 

Chapter IV includes new and original research on whistler wave 
amplification using a wave equation driven by the pendulum equation. As a 
member of a team, | participated in the derivation of the homogeneous and 
inhomogeneous wave equation and electron equations of motion. Next, the 
team developed numerical simulations of VLF wave amplification. Using the 
numerical simulations, the effects of homogeneous and inhomogeneous 
magnetic fields are explored, as well as the effects of a spread in phase velocity. 
The phase- space evolution of the electrons and the growth of the optical pulse 
power and phase are displayed graphically in order to better understand their 
coupled dynamics. 

Much of the research in this thesis has been presented at scientific 
conferences, published, or submitted for publication. The research in Chapter III 
will be submitted to Physical Review Letters under the title "A Quantum 
Limitation on the Start-up of Compact Free Electron Lasers" by K. A. Sturgess 
and W. B. Colson. The research in Chapter IV was presented at the Thirteenth 
International Free Electron Laser Conference in Sante Fe, New Mexico (August 
1991), and resulted in the publication of a paper titled "An Analytical Solution to 
Longitudinal Modes in Free Electron Lasers" by K. A. Sturgess and W. B. 


Colson in Nuclear Instruments and Methods in Physics Research, A318, pages 


576-581, (1992). The research in Chapter V was presented at the Fourteenth 
International Free Electron Laser Conference in Kobe, Japan, (August 1992) 
and accepted for publication in Nuclear Instruments and Methods in Physics 
Peeearch with the title "Short Pulse Evolution in the Boeing/LANL APLE 
Oscillator" by K. Sturgess, D. J. Frost and W. B. Colson. The research in 
Chapter VI was presented at the Fourteenth International Free Electron Laser 
Conference in Kobe, Japan, (August 1992), and accepted for publication in 
Nuclear Instruments and Methods in Physics Research with the title "Analogies 
Between the Free Electron Laser and Whistler Wave Amplification in the 


Magnetosphere” by W. B. Colson, B. Gately, D. L. Caudle, and K. Sturgess. 


ll. FREE ELECTRON LASER THEORY 


A. INTRODUCTION 

An operating Free Electron Laser (FEL) consists of an accelerator capable 
of driving electrons to relativistic energies, an undulator composed of a 
transverse, periodic magnetic field that “wiggles” the electrons, and optical 
devices necessary to collect and transport the amplified radiation. The 
periodically deflected relativistic electrons radiate spontaneously, and in 
addition, couple to the radiation store in an optical resonator leading to coherent 
stimulated emission. The idea of an FEL was conceptually introduced in 1971 
by J.M.J. Madey [2]. Related work which contributed to the FEL include the 
traveling wave tube of Motz [3], and the ubitron developed by Philips [4]. 
Successful experimental operation of the FEL amplifier in 1976 [5] and an FEL 
oscillator in 1977 [6] demonstrated the feasibility and potential of this new laser 
device. 

As mentioned above, the undulator's role in an FEL Is to give the relativistic 
electrons periodic transverse deflections as they travel with the light pulse down 
the length of the undulator. This deflection, provided by either permanent 
magnets or electromagnets, allows the electrons to couple to, and transfer 
energy with the transverse optical pulse along the undulator length, L=WNA, 
where 4, is the undulator wavelength and WN is the number of periods in the 
undulator. An important dimensionless quantity characterizing the undulator is 
the undulator parameter K = eBi,/22mc’, where B is the rms field strength 
over one period of the undulator, e is the electron charge magnitude, m is the 


mass of the electron, and c is the speed of light. In most FEL's, K =0.1 5S. 


There are two main types of undulator, based on the paths electrons take 
within them: helical and linearly polarized. A helical undulator produces helical 
electron paths, which produces and couples to circularly-polarized light. The 
ideal helical magnetic field is constant in magnitude along the longitudinal axis 


of the undulator, and in rectangular coordinates is given by [1] 
B,, = Bl cos(kgz), sin( kz), 0]. (2-1) 


where B is the peak magnetic field strength on this axis, and ky = 27/2 is the 
undulator wave number. The helical undulator design focuses electrons toward 
the undulator axis, since the average magnetic field can be expanded and 
shown to increase off-axis as 


__ Ze 
ae B14 ..), (23) 


\ } 


where ris the distal distance from the axis of the undulator. 

A linearly-polarized undulator produces periodic planar motion in the 
electron beam which in turn produces and couples to linearly-polarized light. 
The linear magnetic field has an rms field strength of B/ J2 which is less than 
the helical case, and so, couples less well to the light and produces half the 


gain. In rectangular coordinates, the linear field is given by 
B, = Bl0, sin( kz), 0} (2-3) 


and provides magnetic focusing in one direction only. 

The accelerators used in various FEL projects range from low energy Van 
De Graff electrostatic accelerators to high energy storage rings, including 
intermediate energy microtrons and linear accelerators. They span an energy 


range giving Lorentz factors y~10—10°, where 77 =,/1-v?/c?, and 


electron beam currents ranging from I*1-+10°A. The typical RF linac 
produces electron micropulses of length, /, ~ mm —+ cm, separated by about 
one to several meters grouped together in a long macropulse of ius to 1ms. 
The transverse radius of the electron beam is on the order of millimeters 
resulting in an electron density of p = 10'* > 10°? cm™?. 

The electron beam quality is characterized by its energy spread and 
emittance. Energy spread describes the spread of electron velocities in the 
beam and is characterized by Ay/y. Emittance is defined as ¢~ 74, where 7, 
is the rms initial position spread of the electrons in the beam, and @ is the rms 
initial angular spread in the beam. Although the angular and position spread of 
the beam can be varied, the emittance of the beam tends to be a constant for an 
accelerator at a given energy. Normalized emittance, ¢, = ye, tends to be 
constant along the accelerator while the electrons are increasing their energies. 

When the electron beam enters the undulator, it is desirable to choose the 
angular and position spreads in the beam so that it does not excessively focus 
or expand as it travels down the length of the undulator. For a fixed value of 
emittance, a small beam radius at one point along the undulator results in a 
large angular spread, and consequently a large beam radius farther along the 
undulator. A small angular spread results in a large beam radius everywhere. 
The condition that minimizes the beam radius along the undulator is called the 
‘matching condition,’ and gives Kkpiy = y@ [1]. This condition determines an 
electron beam radius of r, = Vvé,A/V22K where A=22/k is the optical 
wavelength. A beam matched in this way gives a phase velocity spread of 
Av=2aNKke/y. Since the phase velocity spread must be less than the gain 
band-width, Avs 2z, an upper limit on the allowable emittance is 


e™™ < yA/2nNK. 


There have been many approaches to understanding how the FEL works. 
Among the methods used to develop FEL theory has been a quantum 
mechanical approach, the use of quantum electrodynamics , the application of 
plasma physics theories, and the use of plasma physics many-particle computer 
codes [1]. Early on in the theoretical development, it was shown that a classical 
approach provided an accurate description of the FEL [7]. The classical 
approach using a charged-particle current driving Maxwell's equations will be 


presented here. 


B. ELECTRON DYNAMICS 
The only significant force acting on an electron as it travels the length of the 


undulator is the Lorentz force given by the relativistic equations [1] 


a= ag hE + BX NG, +B, (2-4) 
Cee ee _ 
dt me (2 E,) a 
ae 1 3/80 (6 (2-6) 


where # = V/c is the electron's normalized velocity, B, = B[cos(koz), sin(koz),0] 
is the helical undulator magnetic field, E = E| cos y,— sin y,0| and 
B.=E[siny,cosy,0] are the optical electric and magnetic fields, 
y= kz-at+¢, k= 22/7 is the optical wavenumber, w=ck is the optical 
frequency, and ¢ is the phase of the optical wave. Equations of motion for the 
electrons are obtained by inserting the electric and magnetic fields into the 


equations above. The plane wave form for the electric and magnetic fields is 


accurate near the axis of the undulator. Since there are a total of five equations 
above, and only four unknowns, [7(t),7(t)], one of the equations, the z 
component of equation (2-4), will be ignored. Expanding the remaining 


equations yields 


d vB, — 
_ 1 ale (1- B,)[cos y,— sin y,0] + B,B| - sin{Kpz),COS(koz),0]] (2-7) 





or - = (g,c0s y- A, sin y) (2-8) 
ee: (2-9) 
where £, =[f,,8,,0] is the transverse electron velocity. Since the electrons in 
an FEL are relativistic, v, «c, we have (1—£.)E <<£.B, so that the transverse 
optical force is much smaller than the undulator field force. This removes the 
first term on the right side of equation (2-7). The remaining term can be directly 


integrated yielding 
= —K 
B= = Leos{ ko2), Sin( kez), 0] (2-10) 


The constants of integration are taken to be zero representing an electron beam 
that is injected perfectly into helical orbits. 

The transverse amplitude of the electrons can be found from integrating 
equation (2-10) further. First, assume that y(t) is a constant, which is a good 
approximation for low efficiency FELs. Next, assume that z =£,cf, so that the 
electron oscillation frequency in the undulator is @, =k,fc. The arguments 
within the trigonometric functions of Equation (2-10) can then be rewritten as 


k,z ~@,¢. Integrating now leads to 


X,* Fat [sin wet),cos(eet),0] (2eaah 
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which yields a transverse amplitude of oscillation, Ax, = Ki,/2z7. 
Continuing with pendulum equation derivation, substitute equation (2-10) 


into the energy equation (2-8). Using a trigonometric identity leads to 


dy eEK 
= = sme COS\ (ke + k)z —oat+ g). (2-12) 


The electron's phase is identified in equation (2-12) as €=(k+k,)z-—at. We 
see from equation (2-12) that electron phases such that -—2/2 < (+ ¢< z’2 will 
gain energy from the radiation, while electron phases such that 
nl2 < €+ ¢< 322 will lose energy to the radiation. 

Differentiating the electron phase velocity twice yields e = (k+k,)B.c. 
Inserting the magnitude of equation (2-10) into the definition of y from equation 
(2-9) and differentiating with respect to time yields 7/7 = 772,8,. Combining 


these equations with equation (2-12) gives the pendulum equation 


C= eek ONE cos(¢ + g) (2-13) 
yom 

which governs the electron’s motion in the presence of the combined radiation 

and undulator fields. Defining the dimensionless time as c=ct/L (so that 

t™=0-—>1 as the electrons travel from the beginning to the end of the undulator) 

normalizes equation (2-13) as 


0° 2 
(2) = LE 2ekoKE 


Z 2 


cos(f+ ¢) =|alcos(¢+ ¢) (2-14) 
C Yom 





where () SyeP (oka UE c)*d?(..)/at? and jal = 4aNeEKL/ygmc* is the 
dimensionless optical field strength. The relativistic Lorentz factor y in equation 
(2-13) has been replaced by its initial value 7, in equation (2-14) by assuming 


low efficiency: the energy changes of the electrons are insignificant during the 
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interaction time z. The electron pendulum equation for a helical undulator (2- 
14) is appropriate for low and high gain and correctly follows the development of 
electron bunching. For the linear undulator, |a) should be modified by letting 
K + K(Jy(&)-J,(&)), where J, and J, are Bessel functions and 
E= K?/2(1+K?). 

Using the expressions above, a dimensionless phase velocity can be 


defined by 


X7)= fe LI(k + ko) B, - k]. (2-15) 


The electron experiences resonant undulator and radiation field forces when 
v=0. The resonant relationship between the radiation and wundulator 


wavelengths is 


ee ale’ 


A = “oO Oy (2-1 6) 


Equation (2-16) shows that the operating wavelength of an FEL can be altered 
by changing the electron beam energy ;, the undulator wavelength 4,, or the 


undulator parameter K. 


C. THE SELF-CONSISTENT WAVE EQUATION 

While the Lorentz force law governs the evolution of the electron’s trajectory, 
Maxwell's equations govern the evolution of the optical field in the FEL. 
Spontaneous emission in an FEL oscillator grows to form a classical field with a 
bandwidth proportional to the inverse of the number of periods in the undulator. 
This narrow bandwidth implies that the optical wave varies slowly in space and 
time over one wavelength. This means we may assume (E’ << kE, E << wE, 


¢’<<k¢, ¢<< ad), where (..)'=d(..)/dz, and (..) = d(..)/at. These 


AZ 


assumptions comprise the slowly-varying amplitude and phase approximation 
(SVAP) [1]. 
Neglecting transverse effects, the radiation present in the undulator can be 


taken to be a circularly-polarized plane wave whose vector potential is 
ia ee 
A = sin y,cos v0) (2-17) 


where y=kz—at+¢. Using this vector potential in the Coulomb gauge of the 


wave equation yields 


¥ = aoa 2 “= re ace y,—sin y,0] 





ab ab oa (2-18) 
-2e| + a 2 | sin y,cos y,0| = au 

where the SVAP approximation has been used to eliminate all terms with two 
derivatives. The total beam current density J, is the sum of all the single 
particle currents J, =-ecB,6°(X-7), where 7 is the trajectory of the ith 
electron and B , IS given in equation (2-10). The fast rotating sine and cosine 
functions of w can be eliminated by defining two orthogonal unit vectors 
é,=[cos y,-siny,0| and é, =[siny,cos y,0]. Projecting (2-18) onto these two 


unit vectors separates the wave equation into two first-order scalar differential 





equations 
cE ee eee (2-19) 
eZ tol c 
and 
1 160 2T- . 
Fe aaa ee BON 
[2 | a — 


3 


Making the change in coordinates z=Z+ct in equations (2-19) and (2-20) 
simplifies the differential operator to (f&+c'aa)>c'da. The new 
coordinate z stays with a point on the field envelope traveling at speed c. The 
electrons then drift back with a relative speed o(1 =i Inserting this change, as 
well as the expression for the Current density, yields 


€E _ — 2zeKc 
ak 





cos( w+k)z)d> (X-F) (227) 
and 


Ee _ 27eKC 
5 oie 


Next we substitute y+k,z =¢+¢. 


Silty knack (2-22) 





A volume element much smaller than the coherence volume, but larger than 
one wavelength of light, is selected to sample electrons. Since the fields are 
slowly varying, only one wavelength of phase space is selected and sampled 
from this volume element, because all wavelengths in the volume must evolve in 
the same way [1]. The electrons in one wavelength of phase space are sampled 
to achieve an accurate representation of initial phase space distribution, So the 
sum over all particles in equations (2-21).and (2-22) can be replaced by an 
average over the sampled electrons weighted by the electron density p [1]. 
Indicating the ensemble average of particles at a site z by ) the optical wave 


equations become 





oe = AIK (cos(¢ + 9) (2-23) 
Ct y 
and 
ES - PERCE (sin(c+ $)). (2-24) 


14 


Finally, by introducing the dimensionless time 7zt=ct/L, current density 
j =8N(exKL) p/y2mc?, and field strength |a| = 4aNeEKL/ 2 mc’, equations (2- 
23) and (2-24) further simplify to 


lal=-j(cos(¢+4)),  o= i sinls+ ¢)). (2-25) 
In low efficiency FELs, the energy loss per electron within the undulator is small, 
and so we hold the electron's energy constant at its initial value.of 7. In 


complex notation, the wave equation is 


¢ 


a=- ile“ ), (2-26) 


where a = ale” [1]. Equations (2-25) and (2-26) show that electron phases near 
¢+@ 7 drive the optical amplitude and lead to gain, while the optical phase is 
driven by electron phases near €+¢= 7/2. These equations complete the 
feedback mechanism of the FEL. The pendulum equation describes the 
bunching of electrons which drives the optical amplitude, and the optical 
amplitude controls the amount of bunching that takes place. For the linear 
undulator, the dimensionless current jf is modified by letting 
K -» K(J,(¢)—J,(¢)) where J, and J, are Bessel functions and ¢= K?/2(1+K?). 

The previous derivation of the wave equation disregarded the transverse 
evolution of the fields by making the plane-wave approximation. To include 
these transverse diffraction effects, the three-dimensional FEL wave equation 
uses the SVAP approximation in a similar manner to equations (2-26) and 
results in the parabolic wave equation [8] 


—f= One me 
= vo < jax, = -(je iy (2-27) 
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where V+ =&+&, and transverse coordinates are normalized to the 


y? 
V2 


characteristic mode size, so that x(z/LA)” +x, and y(z/L2)” + y, and the 
longitudinal coordinate is (z—ct)/NA-»z. Equation (2-27) can be used to 


include the effects of electron beam radius and diffraction effects. 


D. FEL DIMENSIONLESS PARAMETERS 

Dimensionless parameters are extremely useful when discussing FEL 
attributes. Several of them have already been defined and used, such as the 
undulator parameter K, the electron phase and phase velocity € and v, the 
current density /, and the optical field a. 

The general method of normalizing consists of finding characteristic lengths 
and times, used as divisors to remove all dimensions. The characteristic time 
scale is given by the time for an electron to traverse the undulator, L/c. The 
characteristic longitudinal length for macroscopic processes is the undulator 
length, L. The characteristic transverse macroscopic length is the optical mode 
waist size [1]. The characteristic microscopic longitudinal distance is the 
slippage distance. or the distance that light passes over an electron bunch as it 
travels down the undulator. Looking at an electron - photon race down the 
undulator, the photon’ beats the electron by a_— distance 
As = Av- At = (c-f,c)L/¢ = (1-f,)N4p. Using (1-,) = (1+K?)/27”, for 
relativistic electrons combined with the resonance condition (2-16), the slippage 
distance can be written as As~ NZ. The slippage distance is the characteristic 
length over which an electron and photon can exchange energy. 

Finally, the Rayleigh length, Z,, is a measure of the optical beam diffraction 
(1] and is determined by the resonator configuration including mirror separation 


and curvature. The optical mode waist is related to the Rayleigh length by 
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mW, =4Z,. The dimensionless Rayleigh length is given by z, = Z,/L [1], and the 
optical mode waist is normalized by ./L4/z so that w, = Jz). 
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lil. QUANTUM LIMITATION ON FREE ELECTRON 
LASER START-UP 


A. INTRODUCTION 

A natural place to begin examining the FEL is at start-up. This chapter looks 
into processes that occur during the first few passes of the electrons through the 
undulator. For special FEL designs, the number of photons in the resonator may 
be few enough to cause quantum mechanical! effects in this otherwise classical 
laser. 

On each pass through the FEL oscillator, the periodic magnetic fields of the 
undulator give the electron a transverse acceleration, causing it to 
spontaneously emit radiation into the resonator modes [9]. On subsequent 
passes through the oscillator, the electrons continue to spontaneously emit 
radiation, but now in the presence of existing radiation. This results in 
stimulated emission and absorption of photons by the electrons. As the number 
of photons in the resonator increases sufficiently to describe a classical field, 
the radiation frequencies with the largest gain grow to dominate the resonator 
mode and the spectrum narrows. If the single-pass gain of the FEL is above the 
resonator losses (the ‘threshold condition"), the radiation power grows 
exponentially to saturation in strong fields. For most operating FELs, the start- 
up process occurs in a few hundred passes. However, a new class of compact 
FELs with short undulator periods 2,, a large number of periods N, and a low 
undulator parameter K = eB2,/2amc’, may experience significant delays in 
Start-up due to large fluctuations caused by the low number of photons 


spontaneously emitted into the resonator mode each pass. In some cases, the 


FEL reaches steady-state without sufficient photons to achieve a classical field, 


and the FEL fails to operate. 


B. COMPACT FEL DESCRIPTION 

Compact FELs are a new, exciting direction in FEL research and offer 
attractive advantages over conventional FELs in several ways. The compact 
FEL has a significantly reduced undulator period, 2) = 1mm, which should 
reduce the size, weight, and cost of the FEL. The shorter undulator period also 
allows production of a short wavelength radiation with a less powerful 
accelerator, since the radiation wavelength is inversely proportional to the 
electron beam energy ymc?, 4 = 4,(1+ K?} eye Unfortunately, decreasing 
the undulator period reduces the coupling of the electrons to the transverse 
radiation field, since the deflection of the electrons in the undulator is given by 
Ax = Ki,/2zy. To make matters worse, there are technological problems with 
obtaining high magnetic field densities in very short undulator periods, thereby 
reducing the undulator parameter to a value K <<7. This problem acts to further 
reduce the gain. The length L=Nz, of the compact undulator must be 
increased to compensate for the reduced gain. 

The accelerator used by the compact FEL to produce short wavelength 
radiation may be much less energetic than that required by a typical FEL since 
the undulator period is reduced. But, because the compact FEL has an 
increased number of periods, the accelerator must have higher beam quality 
consistent with the decreased gain spectrum bandwidth, 1/(2N). In this 
analysis, the beam quality is considered achievable, so that neither emittance 
nor energy spread degrades the FEL gain. In this compact FEL design, the 


undulator period is small which leads to a small undulator parameter K <<7. A 


iS 


large number of undulator periods are used to insure the expected classical gain 


is above the threshold condition. 


C. START-UP OF THE COMPACT FEL 

On each pass through the undulator, the number of spontaneously emitted 
photons from a single electron into the resonator mode is calculated as follows. 
The total power radiated by a single electron is given by [10] 


mezes y 2674 (3-1) 
els 


where B, = Kk,c/y is the electron's transverse acceleration. The energy of a 
photon emitted at the radiation frequency satisfying the resonance condition, 


w@=2y'k,c, is 
E = 2hckyy*. (3-2) 


The number of photons emitted by one electron in one pass through the 
undulator is then obtained by multiplying the power equation (3-1) by the time an 
electron spends in the undulator, At~=L/c, and then dividing by the photon 


energy to obtain 
W, = 2Nak?. (3-3) 


Here a = e*/hc is the fine structure constant. 

The quantity in equation (3-3) is the number of photons emitted into the solid 
angle y~* [10]. As Figure 3-1 shows, this solid angle is typically large enough 
so that most photons miss the mirror at the ends of the resonator and are lost. 


in order to calculate the fraction of the photons captured in the resonator, we 
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must multiply W, by the ratio of the solid angle subtended by the fundamental 


resonator mode. 





Figure 3-1. Diagram showing relative sizes of the emission angle y"' and the largest angle at 
which photons can be Stored within the resonator, 9c. 


The optical mode radius is given by w? =w,’(1+z7/z?), where w, is the optical 
mode waist, 2) = mwe /A% is the Rayleigh length, so that the solid angle 
subtended by the optical mode is ~ 1/ (2Ny”). Multiplying the number of photons 
emitted by an electron per pass, W,, by the ratio of solid angles gives W, = aK ‘ 
[3] as the number of photons emitted per pass that are saved in the resonator 
mode. 

Each electron in the beam emits W, photons into the resonator mode. The 
volume, V, = N*A?>, is used to count electrons and photons and is determined 
by a longitudinal length corresponding to the gain bandwidth NA, and the 
transverse mode area Li =NZ2/, [1]. The number of photons spontaneously 
emitted in one pass into this volume is S. = pFak*N’#2, where p is the electron 
beam particle density, F = zr,” / mw, is the filling factor describing the transverse 
overlap of the electron and optical beams, r, is the electron beam radius, and 


Wo =(LA/3n)"" is the radiation mode waist . 


Z| 


After a few passes, there are N, photons in the volume V,. The resulting 
radiation field follows a Poisson distribution giving a fluctuation in the photon 
number of N°”? [11]. If these fluctuations are larger than the gain bandwidth of 
the F EL, 1/(2N), they act to reduce gain. To include the effect of these photon 
fluctuations we consider a simple analogy. Since the dimensionless radiation 
field is given by la! = 42NeKLE/y5mc?, where E is the radiation’s electric field, a 
fractional fluctuation in E is equivalent to a fractional fluctuation in the undulator 
parameter K. Thus, to the electron, the radiation fluctuations are equivalent to 
an imperfect undulator field of Ni’? coherent periods. Clearly, N,”* = N 
signals the presence of a classical field within the resonator. The gain in the 
undulator is proportional to the dimensionless current density 
— BN? (e2Ki9)° pF / yamc* , although j must be corrected by the ratio 
(Ni/2 /n) for the calculated gain to account only for the coherent number of 
periods in the undulator. However, there are (N/N,”?} coherent undulators 
acting in series. Thus, the actual gain per pass after one pass, including 
quantum fluctuations, is estimated as 
N, 


Ge Co 7 2 


(3-4) 


where Gp = 0.135,/ is the classical small-signal gain. If the number of photons 
spontaneously emitted on the first few passes is low, the photon fluctuations may 
be large and lead to a significant decrease in gain. 

After many passes, the number of photons in the resonator will change due 
to continued spontaneous emission, the resonator Q losses, and possible gain 
from stimulated emission. The differential equation describing the evolution of 


the photon number, N_, over many passes, N, is [11] 


Ze 


dN 
= = S, + N,(6-Q"), (3-5) 





where Q represents the loss within the resonator. The FEL may remain below 


threshold, G<Q™', and yield a steady state when the number of photons 


N2 48.6.0? / \"* 
N. = Vee ||| eweee ey, 
Pa IeG N2} |. (3-6) 


Equation (3-6) shows the possibility of a complex number of photons, which can 


reaches 





be interpreted as the absence of a steady-state solution and no start-up of the 
FEL. From this interpretation, it is inferred that the start-up condition is 
4S,G.Q?/N? >1. Substituting S, = pFaK?N?772,, the dimensionless current 
density /, the small signal gain G, =0.135/, and introducing the Compton 


wavelength 7, = h/mc leads to a general start-up condition 
0.135( jQ)* ay > (22)° NX. (7) 


In the classical limit, A, >0, all FELs satisfy this requirement and start. A 
typical FEL (such as the original Stanford SCA/FEL designed with j/Q =(4z)’, 
Lorentz factor y=80, 2=3um, and N=160) satisfies the requirement of 
equation (3-7) by two orders of magnitude. 

On the other hand, a compact FEL with parameters of y =24, A= 86um, 
jQ=70, and N=10° results in no start-up. Two figures, 3-2 and 3-3, appear 
below. Figure 3-2 represents the results of a conventional (‘'classical”) 
calculation which does not include quantum fluctuations. This calculation and 
that figure erroneously predict successful start-up. Figure 2-3, which does 


include quantum fluctuations, correctly predicts that the FEL does not start up. 
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The numerical simulation including the classical Spontaneous emission line 
spectrum is used to determine the photon evolution per pass without the 
quantum fluctuation term included in the spectrum. The bottom panel shows the 
classical spontaneous emission spectrum while the middie panel shows the 
photon number evolution over n=300 passes with the final spectrum at the top. 
Over a few passes, the photon number spectrum shifts to the maximum gain 
frequency and narrows to a classical FEL line spectrum. The gain spectrum Gy 
is displayed with the dotted line showing the loss per pass. The total number of 
photons spontaneously emitted per pass is given by S, =278. The total number 
of photons in the final laser spectrum is N,, = 2.9 x 10'°. Normally, high-power 
Saturation would limit the final number of photons to a value less than N,, before 


n =300 passes, but high-power saturation is not included in this simulation. 


on (Vv) =S,,(V) +N, (V) [G (Vv) -1/9] 









Coherence Evolution 
j= 0.813 Q=90 
K=0.01 N=1000 






S.=278 


G= 0.11 
Y 






Figure 3-2. Photon number evolution using small-signal gain. : 
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Figure 3-3 displays the results of including quantum fluctuations in the gain 
term. After n=2000 passes, the upper panel shows the photon number spectrum 
begins to increase in amplitude and evolve toward the maximum gain frequency. 
Before reaching the threshold, a steady state is reached with N,, = 20: 
photons, where the increase in photons per pass due to spontaneous emission 
S, =278, and stimulated emission GN, is equal to the losses in the system, 
QN,. The actual gain of the laser, including the quantum fluctuation term, Is 
below the losses of the laser, as shown in the center panel. Although this 
compact FEL would be expected to work according to classical FEL theory, 
Figure 3-3 shows a steady state photon number is reached in the resonator 
below that required for a classical field’) Thus, the gain never exceeds the 


losses and the laser does not start. 


ON, (V)=S.(V) +N, (Vv) [G (Vv) -1/Q] 
G(v) = 0.1355 (N, /n*) 








Coherence Evolution 
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Figure 3-3. Photon number evolution using the modified equation (3-4) as the gain. 
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D. ANALYTICAL RESULTS 

Assuming a form for the gain of G = GN, liad the photon rate equation (3- 
5) can be solved analytically for the photon number evolution per pass. Again, 
this solution is valid only up to N, = N* at which time the model used to develop 
the gain formula is no longer valid. The solution separates into two cases: one 
when the start-up condition (3-7) is satisfied, and the other when it is not. If the 
start-up condition is not satisfied, the photon number after the n” pass is given 


by 
Ne Nt 1+ {1 + a (3-8) 
” — 2G,Q 65 Soa 


where q = (1- 4G,S,Q@?/n?)* and r=|(1+q)/(q-1)}. This solution displays a 
linear growth dependence on n for small n. However, there are never enough 
photons in the resonator for the gain to overcome FEL loss defined by Q, soa 
steady-state is reached with an insufficient number of photons in the resonator 
to form a eee field. If the start-up condition (3-7) is satisfied, the photon 


number after then” pass is 





2G,Q 1+ 2 tan{ 4 , (3-9) 


where p = (4G,S,Q? /N? - 1). This solution for the photon number evolution 
grows slowly per pass initially due to contributions primarily from spontaneous 
emission. As the number of photons in the resonator grows, stimulated emission 


begins to add photons and their number increases dramatically. This continues 
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until the fluctuations in photon number are much less than the gain spectrum 
bandwidth resulting in classical gain. 

The threshold for which stimulated emission produces as many photons as 
does spontaneous emission is obtained by setting S, =N,G. Dattoli [12] gives 
an expression for this threshold number of photons as (W,). = ya] (3277Ni . } 
using a quantum mechanical formalism. This expression can be reproduced 
with all the correct functional proportionalities much more simply by setting the 
gain equal to the small signal gain G=0.135/. However, when quantum 
fluctuations are included, and the gain is determined by (3-4), the condition for 


the stimulated emission to equal the spontaneous emission is 





N32 V/2 
,),=( ae] @-10) 


This new condition is different from Dattoli's and takes into account the reduced 
gain initially present due to fluctuations in the radiation field. It can be seen that 
if the laser start-up condition is not fulfilled, the steady-state photon number 
within the resonator, N, will be less than the threshold number of photons 


required for equal spontaneous and stimulated emission rates. 


E. CONCLUSIONS AND FURTHER RESEARCH SUGGESTIONS 
The type of FEL that may exhibit start-up problems due to quantum 
fluctuations is one that attempts to produce short optical wavelengths with a low 
energy electron beam. The low energy electron beam necessitates a small 
period undulator, which in turn reduces the undulator parameter. Both of these 


properties act to reduce the expected gain of the FEL. In order to achieve 


2/ 


sufficient gain, the number of periods in the undulator must be increased. 
However, this design logic can lead to an FEL that may not operate due to 


quantum fluctuations during start-up. 
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IV. ANALYTICAL DESCRIPTION OF LONGITUDINAL 
MODES 


A. INTRODUCTION 

The optical power in the FEL oscillator starts from spontaneous emission 
and grows to strong-field saturation over many passes. During the crucial 
Startup process, coherence develops in the weak optical field and small 
microscopic bunching develops in the electron beam. In this chapter the FELs 
considered are typical in the sense that the quantum effects previously explored 
do not occur. The FELs considered here have weak classical fields, low gain, 
and electron pulses much longer than the slippage distance NZ. 

Assuming low gain, a long pulse, and weak optical fields, Elleaume [12] 
showed that FEL optical pulse evolution could be described by a Schroedinger- 
type equation with a non-Hermitian harmonic oscillator Hamiltonian . Using the 
same assumptions, Dattoli [13.14,15,16] employed methods of quantum optics and 
an integro-differential equation to obtain his solution to the same problem. In 
each case, the mathematical result involves second-order derivatives. In this 
chapter, we expand the optical and electron pulse shapes in a Taylor series to 
derive a first-order partial-differential equation describing the optical field 
evolution over many passes through the undulator. This equation is used to 
investigate coherence development and the effects of the electron pulse shape 


on the optical field amplitude and phase. 


Zo 


B. MATHEMATICAL DEVELOPMENT 

In order to generalize the optical field in the longitudinal direction, to allow 
for multiple optical modes, we must follow the optical field envelope evolution all 
along z. The extension to multi-mode theory by transforming a-—> a(z) is 
equivalent to allowing the optical field a continuum of wavenumbers a -—> a(k). 
Fully extending the theory uses a- a(z), v-> Kz), and €> d(z).[1] As 
discussed in Chapter ll, all longitudinal distances are normalized to the slippage 
distance, the distance over which the light and electrons can exchange 
information, so that z/NA — z. The electrons slip back one slippage distance, 
or Az = 1, with respect to the light during the interaction time z=0 —7 leading to 
NZ = Aza 


The FEL interaction can be described by the optical wave equation [1], 


a(z,) = jo(z + 2)(exp[-i¢(z + z,7)]) (4-1) 


coupled to the electron pendulum equation [1], 


dz, t) = Az, t) = la(z — f, t)\cos| ¢(z, tr) + @(z— 7, )|, (4-2) 


where oy =) ce ae k))z- at is the’ electron '§ phase, 
v=L[(k+k,)f, —k] = ¢ is the electron phase velocity, la(z, r) is the optical 
electric field amplitude, dz,t) is the optical phase, and <..> represents an 
average over all of the electrons in the beam. The electrons’ pulse shape is 
described by o{z) with peak value of unity. The longitudinal field sites z refer to 
positions on the optical field envelope. In (4-1), the light remains fixed in z while 
the electrons slip back to a field site z + 7, while in (4-2), the electrons remain 


fixed and the light slips ahead by z — z as it moves through the undulator. This 
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slippage allows the light and electrons to exchange information, and allows the 
electrons to transfer of information from one optical site to another.[1] 

When the FEL is operating in the low-gain regime, we expand equation (4-1) 
to lowest order in the current density, /, and assume the field changes are small 
over a single pass, al <<{al and o<< g. Since the optical fields are weak, the 
electron phase may be expanded to lowest order in the field magnitude lal, 
giving ¢=¢, + r+... . Here Z, is the initial electron phase, and v, is the initial 
electron phase velocity. The long pulse shapes may be expanded in a Taylor 
series in the normalized slippage distance,r: o{z+7)=o0(z)+ 70'(z)+..., 
la(z- z)|=|a(z)j- qa(z)/+... and  ¢(z-7) = ¢(z)- 7¢'(z)+... where 
(..)'=d(..)/dz. These assumptions allow the expansion and integration of 
equation (4-2), 

(z+72,t) = f+ %,+ ly cos(¢, + Vott g) + COS(6y + ¢)- vorsin(gy + ¢)| 


Vo 


s 


; a 3 |2sin( Zo + vot+ $) - 2sin(o + ¢) -—2vyrC0S(Cy + $) 9 7 (% + ¢)| 








Ae 

+ ag /2 cos(¢o + vot + ¢) -— 2C08(cq + $) - 2vorSin( fo + $) + Vo 7 cos(Ly + ¢)] 
0 

Foeeeeey aie 


where |al=|a(z)|, lal =|a(z)|, ¢= o(z), and ¢" = ¢'(z). 

For a monoenergetic beam (where all electrons have the same initial phase 
velocity, v,) the aeen dynamics expressed in equation (4-3) can be 
incorporated into the optical field equation (4-1) using the phase average 


2a 
ei I, (..)déy / 22. The complex optical field equation then becomes 
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da | ei + fame tas) ie a 
Or Vp v4 2 Vp (4-4) 


he en!%r bn iVo — 


Mo 


viaje 
where a =|a(z)\exp(id{z)). This equation includes the electron dynamics of 
equation (4-2) using the assumption of low gain, weak fields, and long pulses as 
expressed in equation (4-3). 

Integrating over a single pass, r=O-—7, yields the following partial- 


differential equation for the net change, Aa, in the complex optical field, 


Aa = jacg( v)- = ja =H v)- 5 jo=9'( v9) (4-5) 


with constant coefficients for a given initial electron phase velocity 


2-2e -iy(t+e}  6-6e —Gine* + 4(14 267} dg 
g = ee? hh = se? 
_ Dia 2M Avy 


With no pulses, o=7 and a’=O so there are no longitudinal modes. The 
second and third terms on the right side of equation (4-5) are zero. In that case, 
equation (4-5) describes the usual single-mode, low-gain evolution of the 
complex optical field [1]. The real part of g(v,) is the familiar antisymmetric gain 
spectrum in v,, and the imaginary part describes the symmetric shift in the 
optical phase [1]. The phase velocity v, determines the resonance between the 
optical wavelength, undulator wavelength, and electron z velocity; maximum gain 
occurs for v, = 2.60616. The functions h(v,) and g’(v,) in (4-5) determine the 
effects of slippage and the electron pulse shape on the optical field evolution. 


Elleaume [13] and Dattoli [14,15,16,17] obtained similar functions of vy, in their 
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treatment of the long pulse problem, but they also keep some second-order 
derivatives in the expansion of the field. 

Over many passes n through the undulator, the effects of resonator losses 
and desynchronism must be included [1]. Desynchronism measures the 
mismatch at the beginning of the undulator between the arrival times of the 
succession of electron pulses from the accelerator, and the rebounding optical 
pulse in the resonator cavity. The dimensionless desynchronism is 
d=-2A4S/NA, where AS is the difference between the actual resonator length 
and the synchronous cavity length [1]. The resonator loss over many passes is 
described by the resonator quality factor Q. The resulting first-order partial 
differential equation gives the development of longitudinal modes in a long pulse 
over many passes n in weak fields of a low gain FEL, 

ca 0a a Ca 


ap = jaog( v9) - 5 ja = h( vp) - Pees Coa 9 oe (4-7) 


The Peaplicated optical pulse evolution determined by equation (4-7) depends 
only on the dimensionless current /, the resonator Q, the desynchronism d, the 
electron pulse shape o(z), the initial optical pulse amplitude |, (Z) and phase 
¢,(z), and the initial phase velocity 1. 

Since the shape of the phase, 4(z), is free to evolve, the phase velocity y is 
somewhat arbitrary. For instance, if the phase has a linear slope, gz) = Akz, it 
iS equivalent to a change in wave number, k >k +k, and alters the effective 
value of the phase velocity. Therefore, the phase velocity v, in equation (4-7) 
can be chosen for convenience with no loss of generality. A unique and 
important value of y is the phase velocity giving peak FEL gain, g.'(¥) = 


The optical pulse will naturally evolve towards v5 =2.60616 because of mode 
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competition. In fact, vy, cannot be chosen far from v,, because the slope of the 
phase required to evolve to peak gain can exceed the restrictions of the slowly- 
varying amplitude and phase approximation. This phenomenon can be seen by 
beginning with the phase velocity v, = L[(k+k,)@, - k] and deriving the change 
in phase velocity for a change in wave number Av=-Ak=-g¢’. It can be seen 
that the slope of the phase provides the change in wave number. Next, using 
the slowly-varying phase approximation, the change in the slope of the phase 
must be small as the optical field makes one pass in the undulator, so AvAr< 1. 
During one pass, the optical phase interacts over a dimensionless distance of 
Az which can be linked to the size of the electron pulse o,. This in turn gives a 
new requirement AvAo, < 1. Since one of the basic assumptions of this analysis 
that o, >> 1, the restriction Av<<1 is implied. When vy, =v, and g,'=0, the 
equations for the amplitude and phase in (4-7) become uncoupled, and are 


given by 


C a! 2 a , a 
ae ee 4 igo) =(jg,o+4jhio'-Y2Q)ial, (4-8) 


6b feta wae | 
= t(d- $i) = = jgo- 3 iho’. (4-9) 


The equations describing the evolution of the optical amplitude and phase can 


now be solved separately. 


C. ELECTRON PULSE EFFECTS | 
In this section, three different shapes for the electron pulse are considered 
and then used to solve equations (4-8) and (4-9) analytically. The method of 


characteristics [18] will be used to solve equations (4-8) and (4-9). 
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The left-hand-side of equations (4-8) and (4-9) can be seen to be the 
directional derivative of ja(z,n)| and ¢{z,n) respectively, in the direction of the 


vector |1,(d - 4 jg’c)|. The equations 


=~ 2 (4-10) 
Os 2 
and 
dn _ 
4-11 
es ( ) 


determine a family of curves z=2z(s) and n=n(s) whose tangent vector, 
[z’(s),n’(s)], coincides with the vector |1,(d-4jg’c)]. Therefore, the derivative 


of a(z,n) along these curves becomes 


oja| _ dla[z(s).nls)] _ Ala az | 2 lal en _ 
ds Os - 6z @S én és 


4-12 
é lal. +(d- ighe is) “8 = (i, + iho’ | le a 
én 2) a2 > en 
and 
dg _Agz(s).n(s)] ap ez | ap én _ 
ds ds Oz €S oOnds (4-13) 


A iho’ 
#4 fa- igio ie \ % ob -( ia, + jh. fg) 
On 2 te 2 2 Se 


after using the chain rule and equations (4-10) and (4-11). The family of curves 
z=2(s), n=n(s), and a(s) = ) = |a(s s)ie'**, determined from the solutions of the 
system of ordinary differential equations (4-10) through (4-12) or (4-10),(4-11), 
and (4-13) are called the "characteristic curves" of the partial differential 
equations (4-8) or (4-9) [18]. The approach of solving equations (4-8) or (4-9) by 
using the solutions of potetene (4-10) through (4-12) or (4-10),(4-11), and (4- 


13) is called the “method of characteristics". The existence and uniqueness 
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theorem assumes certain smoothness conditions on the functions of z and n 
which guarantee that exactly one solution curve (z(s),n(s),a(s)) of equations (4- 
10) through (4-13) passes through the initial condition point (z,,M),a)). Using 
the initial values to create the parametric initial solution curve gives z=z,, n=0, 
and al = la, (z,) where |a,(Z,)| is the initial optical field shape and ¢= ¢,(z,) 
where ¢,(z,) is the initial optical phase. Equations (4-10) and (4-11) can be 
solved for z and rn as functions of s and z,. These solutions are then inverted 
and equations (4-12) and (4-13) are solved and substitutions made to obtain 
la(z,n)| and ¢{z,n). The three specific cases follow. 
1. Continuous Electron Beam 

For the special case of a continuous electron beam, o{z) =1 for all z, at 
the phase velocity for the peak gain, vy, the optical pulse can still evolve over 
many passes. The first-order, partial-differential equations (4-8) and (4-9) are 
solved using the method of characteristics. 


Solving equations (4-10) through (4-13) yields 


2(s,Z9) = Zz +(d - 3 jgi)s. (4-14) 

n(s,Z)) =S, (4-15) 

a(s,z0) = [a(zoexe| (i 3} 5| aul 

and WAS,Zp) = 9S + ¢o(Zo) (4-17) 


which can be inverted to solve for s and z, in terms of n and z. These results 
are then substituted into the expression for jal and ¢ to yield solutions for the 


optical amplitude and phase, 


la(z,n)| = |a (Zo )/exp[ (ig, - /2Q)n] , (4-18) 
and | 
G(Z,n) = j9in+ do(Zp) , (4-19) 
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where |a,(z,)| and ¢,(z,) describe the initial optical pulse, g, = Re(g) is the real 
part of g(%), g,=Im(g) is the imaginary part of 9/ vo), and the position 
Zz, =z-(d-4/g,')n describes the optical pulse centroid shifting over n passes 
through the resonator due to desynchronism. lf d <2 jo. . the optical pulse shifts 
back over many passes toward the trailing edge of the electron pulse. Further, if 
d>i jo, . the optical pulse shifts forward toward the front edge of the electron 
pulse. The shape of the optical pulse amplitude and phase remains the same as 
the initial pulse throughout the evolution. Figure 4-1 shows a SSisa of the 
analytical solution (4-18) and (4-19) with the results of a numerical simulation for 
the same conditions. The electron pulse is continuous in length with all 


electrons at phase velocity vy for maximum gain. 
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Figure 4-1. Optical pulse evolution comparing a numerical simulation with the analytical solution 
for a square electron pulse. 
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The desynchronism is d=0.018; the initially weak optical field amplitude is 
parabolic with a maximum of a, = 0.001 and a full width of o, =10; the initial 
phase is ¢,(z)=0 for all z; and the resonator loss is determined by Q=100. 
The lower-left panel displays the continuous electron pulse shape with o=1 for 
all z. The middle-left panel shows the evolution of the optical pulse over n= 60 
passes with the final optical field amplitude and shape in the top-left panel. The 
analytical solution gives a slightly larger amplitude than the numerical 
simulation, but the centroids and shapes of the two solutions are in close 
agreement. The middle-right panel shows the evolution of the optical phase 
over the n=60 passes with the analytical and numerical solutions after the final 
pass displayed at the top-right. The similarity of the two solutions is quite 
apparent within the portion of the window occupied by the optical pulse with 
significant amplitude. In regions where the optical amplitude is small, the optical 
phase becomes meaningless. 
2. Parabolic Electron Beam 
The parabolic electron pulse has the form o{z) =1-2z*/ o,* for 
z| = Ce) 2, and o{z)=0 otherwise. The parabolic electron pulse shape is 
often used in FEL simulations. In addition, the parabolic equation expresses the 
Same z dependence as the first two terms in a Taylor expansion of any pulse 
shape. Using the parabolic electron pulse shape, equations (4-8) and (4-9) can 
be solved analytically for arbitrary initial optical pulse amplitude la (Z)|, and 
phase, ¢,(z) using the method of characteristics. 
The solution to equations (4-8) and (4-9) yields results for two cases 
depending on whether the optical pulse shifts forward due to large 


desynchronism, d >3/9,;', or backward due to large gain, d<iyg,’. 
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In the first case, where desynchronism dominates gain, d >4/g,’ and the 


optical pulse tends to shift forward, equations (4-10) through (4-13) yield 


s=— ton #520) [)ton-( 22 [8)) (4-20) 
| big’ OC, b o,\V D 


Alcea) = s. (4-21) 








la(s, Zp )| = a:(z)en| | = 55 i 2a 5 
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a cos(q) 
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+ 
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(4-23) 
2h 


where b =d —4/jg,' and 


eee 
g = tan” ign z* BE e(big hs. 
o7b cai 


Substituting equations (4-20) and (4-21) into equations (4-22) and (4-23) gives 
the optical amplitude 
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a(2z,)| = es(zaler (1, f+ 202) | 


20 gy 
' v2 \1 (7979 
cs (25) | 
oO 
ex 8 -29,2 + 218 tan(q) eee. Se ae (4-24) 
Gi (jg, N72 cos(q) : 


and the optical phase 
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(4-25) 
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1/2 
where Zp) = (o,7b / 4G; ‘| tan(q), and q is redefined as 


5 1/2 
q = fan (42 |-Ztve )' . 
Geno op 


In the second case, where gain dominates desynchronism, d <+/Jg,’ and 


the optical pulse tends to shift backward, equations (4-10) through (4-13) yield 


s.20) (Bi ao 
§ = — 22 fp.| ee | ee | 
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(4-28) 
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where w = 2s((—b)(jg,'))”* / oy. 
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Substituting equations (4-26) and (4-27) into equations (4-28) and (4-29) 
and redefining w = 2n((-b)(jg;"))” / a, yields the optical amplitude 








1/2 
| File? Def ee 
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and the optical phase 
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The upper sign is used when |z| <—bo, / (—bjg,’)”*, and the lower sign is used 
when z is not between these limits but is within the limits of the electron pulse. 
Outside the limits of the electron pulse, the phase remains constant and the 


optical amplitude has the form 
la(z,n)| = (1 -2Az- dn)? /o?) exp(—n/2Q) (4-32) 


where oa, is the width of the optical pulse. 

The expressions (4-24), (4-25), and (4-30) through (4-32), provide an 
analytical solution for an initially parabolic optical and electron pulse. Figure 4-2 
shows the comparison of the analytical solution with the results of a numerical 
solution for the same conditions. The desynchronism is d=0.015, the initial 
field amplitude a, = 0.001 for all z, the phase ¢,(z)=0 for all z, the electron 
pulse length is o, =10 at phase velocity vy for maximum gain, and Q=100. 
The initially long optical pulse takes on the apparently Gaussian shape shown at 


the upper left after n= 60 passes. 
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FEL Pulse Evolution **** 
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Figure 4-2. Optical pulse evolution comparing a numerical simulation with the analytical 
solutions for a parabolic electron pulse. 


The analytical solution gives a slightly greater amplitude than the numerical 
simulation after n=60 passes, but the shape and centroids of the pulses are 
nearly the same. The difference in amplitude is probably due to an 
accumulation of small numerical errors in the calculation of gain each pass. A 
difference in gains of only 0.3% over n=60 passes would account for the 
discrepancy. 

When d=0, the second case is appropriate because g,'( 4) > 0 and 
equation (4-30) shows that the field ja(z,n)| > 0 as n>. This situation is the 
well-known result that a pulsed FEL cannot operate at zero desynchronism [1]. 
Physically, the centroid of the optical pulse is shifting back over the electron 
pulse by an amount 3/9,’ for each pass through the undulator. Eventually there 


is no overlap between the pulses and the optical field decays to zero. Figure 4-3 
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shows a gain versus desynchronism curve for the parabolic electron pulse and 
initially constant optical pulse. The desynchronism for maximum gain is located 
at d =4/J9,', which is the desynchronism necessary for the optical pulse centroid 
and electron pulse centroid alignment at the beginning of the undulator each 


pass. 


a2 
Desynchronism 
Figure 4-3. Gain vs. desynchronism for a parabolic electron pulse. 





3. Triangular Electron Pulse 

Another electron pulse shape that can be solved analytically using the 
method of characteristics with equations (4-8) and (4-9) is the triangular pulse 
shape of the form ‘o{(z) =1-|z|/o, with width o, for |z|<o, and o{z)=0 
otherwise. This case compares the differences in optical pulse shape and 
phase caused by the z* dependence of the parabolic electron beam with the z 
dependence of the triangular electron beam. Due to the discontinuity in the 
triangular pulse at its peak, z =O, the solution has slight differences in form for 


z <0 and z->0O. Anerror near z = 0 is to be expected, since the triangular pulse 
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shape violates the assumption that the electron pulse shape varies slowly; 
o’(z) <<1. As inthe parabolic case, the optical pulse slipping backward due to 
large desynchronism, d > jg,;', or forward due to predominant gain, d >4yq,’, 
appears. However, the solutions to equations (4-8) and (4-9) for the triangular 
Breeton pulse shape do not lend themselves to being displayed in a form that 
makes these two cases obvious. Therefore, we shall solve for cases of 
Iz| < 2fo, / jg;' and Iz| > 2fo, / jg;’. 

Solving equations (4-10) through (4-13) for the first case of |z| < 2fo, / jg;' 
yields 








z(s,Z)) = = _f + Mi*o os} 2) i (4-33) 
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where f =4/g,'-d. Substituting equations (4-33) and (4-34) into equations (4- 
35) and (4-36) gives the optical amplitude 
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The second set of analytical solutions are for |z| > 2fo, / jg;'. Solving 
equations (4-10) through (4-13) for the case of |z|>2fo, / jg,’ gives 
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Substituting equations (4-39) and (4-40) into (4-41) and (4-42) gives the optical 


amplitude 


la(z,n)| = as(zaleno (28° — SL _ 5) 
{ Zz 


x exp— (423 
I9; 


Bes IGjZ 
2 
where |a,(z,)| is the arbitrary initial optical pulse amplitude, and the optical 
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where ¢,(z,) is the initial optical phase and 
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For both the amplitude and the phase solutions, the upper sign is used when 


—_f+ IG; z 
Oe 








z >0O and the lower when z <0. 

The expressions (4-37), (4-38), (4-43) and (4-44) provide an analytical 
solution to the optical pulse amplitude, shape, and phase in the presence of a 
triangular electron pulse. Figure 4-4 shows a comparison of a numerical 
simulation and these analytical solutions for the same initial conditions. The 
desynchronism is d = O, the initial field amplitude is a, =0.001 for all z, the 
phase is ¢,(Z) = 0 for all z, the electron pulse length is o, =10 at phase velocity 
v, for maximum gain, and Q@=100. The lower-left panel shows the triangular 
electron pulse at the beginning of the undulator wnen z =O, and at the end when 


t=1, showing the slippage distance. The middle-left panel shows the optical 
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pulse evolution over n= 40 passes with the final pass displayed at the top-left. 
The analytical results are slightly larger in amplitude and more sharply peaked 
than the numerical simulation, but the centroid and shape correspond well. The 
middle-right panel shows the optical phase evolution over n= 40 passes with 
the final phase displayed at the top-right. As can be seen, there is excellent 


agreement between the analytical and numerical solutions. 





Figure 4-4. Optical pulse evolution comparing a numerical simulation with the analytical 
solutions for a tnangular electron pulse. 


D. COHERENCE DEVELOPMENT 
The evolution of the optical amplitude and phase over many passes can be 


used to analytically examine coherence development. Taking the parabolic 
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electron beam with d <4 jg' as an example, the optical amplitude solution (4-30) 
on the first pass is given by |a(z,n)| =a (Zo)| where |ao(Z>)| represents an 
arbitrary optical pulse shape as a function of z)(z). After many passes,n— 0, 


the expression for Z) becomes 


OOF 








“Oa n> (4-45) 
(—big, | 
and the optical amplitude becomes 
a\Zz,n =lalZ ex r =e 25h 9 j ” 
| ( ) | ( 0) p 9 20 J [28-8 
‘gs (4-46) 
x@Xp) + -49,0,| 42.52] 4 r(z) + 4” 
19, Fiza) ea | 


To obtain (4-46) from (4-30) we note that each term involving w as a negative 
exponent goes to zero since wo as no. Since Zp, is a constant after 
many passes, the first term in (4-46), which represents the effect of the initial 
optical pulse shape on the pulse shape after n-> o, has been reduced to a 
constant. This indicates that, after many passes, the optical pulse shape is not 
affected by the characteristics of the initial optical pulse. Therefore, this theory 
predicts that an initially random optical pulse evolves into the same shape as an 
initially smooth optical pulse. It follows that the theory should be applicable to 
FELs starting from classical spontaneous emission through coherence 
development to the onset of strong optical fields. It should be noted that the 
spontaneous, random field described here is classical, and does not involve the 


quantum processes analyzed in the previous section. 
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To observe the evolution of coherence development in a more quantitative 
way, several numerical solutions of equation (4-7) were obtained. Shown in 
Figure 4-5 is the coherence development of an initially random continuous wave 
optical field with average strength a, =1, dimensionless current /=1.8, initial 
phase velocity 4, and parabolic electron pulse shape o{z) =1-22z7/ o,” for 
\z|< V2, and zero for |z|> /20,, of width o, with losses determined by Q =6. 
The desynchronism is d=0.02, and the initial optical amplitude has a random 
component added at each size z. The plot of j(z—7z) at the lower-left panel 
shows the parabolic electron pulse shape and the small slippage as rt goes from 
zero to one. The optical amplitude evolution la(z,n) in the middle-left panel 
shows an initially random pulse which becomes much more coherent over n= 40 
passes. The final pulse shape in the upper-left panel is much smoother and 
shows a high degree of coherence. The gain spectrum, g,(v) in the lower- 
center panel shows the antisymmetric gain function for reference. The evolution 
of the logarithm of the optical power spectrum, log (P( v,n)) is shown in the 
middle-center panel. The evolution of coherence leads to a spectrum centered at 
the phase velocity for maximum gain. The final power spectrum in the upper- 
center panel is narrow indicating long range phase coherence. The optical 
phase evolution ¢{z,n) in the middle-right panel shows the development of the 
optical phase over n= 40 passes with the final phase profile at the upper-right. 
The power evolution P(r) in the lower-right panel shows power initially decrease 
as incoherent parts of the optical field are removed by losses. Then the power 
begins to grow as the pulse becomes more coherent and more centered at the 
frequency for maximum gain. Figure 4-5 also implies that the optical shape is 


determined primarily by the shape of the electron pulse. 


SO 


** Analytic FEL Pulse Evolution ** 
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Figure 4-5. Optical coherence development from a random pulse solving (4-7) numerically. 


E. CONCLUSIONS AND FURTHER RESEARCH SUGGESTIONS 
By assuming a long electron pulse, weak optical fields, and low gain, the 
electron and optical pulse shapes have been expanded in a Taylor series in 
order to solve the pendulum and wave equations. The result is a single first- 
order, partial-differential equation describing the evolution of an optical pulse. 
By including the effects of desynchronism and resonator losses, a more 
complete partial-differential describes optical evolution over many passes in an 
FEL oscillator. As a direct consequence of deriving the differential equation, the 
familiar antisymmetric gain.curve and symmetric optical phase shift are obtained, 
as well as functions describing the effects of slippage and the electron pulse 


shape. Solving the differential equation (4-7) for a continuous electron beam 
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yields the simple single-mode gain value of 0.135/ and results in no change in 
the optical pulse shape. Electron pulses with parabolic and triangular shapes 
are shown to yield gains less than predicted by single-mode gain calculations 
since there are fewer electrons distributed over the entire length of the electron 
pulse. The shape of the electron pulse is shown to determine of the shape of 
the optical pulse in steady state after many passes. Coherence evolution of an 


initially random pulse is also studied. 


S2 


V. SHORT PULSE EVOLUTION IN THE BOEING 
APLE OSCILLATOR 


A. INTRODUCTION 

The previous chapters have looked at FELs with low gain and weak optical 
fields. We now examine an FEL oscillator that can produce high gain and strong 
fields, but is designed to produce only moderate optical fields as a seed laser 
for an FEL amplifier. Indeed, the problem investigated is how to make an 
oscillator work poorly, or to reduce the electron-light coupling, so that only a 
small energy spread is induced in the electron beam. This chapter discusses 
the Boeing/LANL Average Power Laser Experiment (APLE) which is designed to 
produce 100 kW of average power for 3 minutes from a FEL using a Single 
Accelerator, Master Oscillator, Power Amplifier (SAMOPA) design, shown in 


Figure 5-1 [19,20]. 


1. Photo Injector 5. Bunchers 
2. Accelerator Section 6. Accelerator Section 


3. FEL Oscillator 7. FEL Amplifier 
4. Accelerator Section 





Figure 5-1. Schematic diagram of the Boeing APLE SAMOPA 


33 


In Figure 5-1, the accelerator section (2) preceding the FEL oscillator (3) 
gives the beam 1/7 MeV energy with an instantaneous energy spread of 0.3% at 
an average current of 0.23 A and a normalized emittance (90% edge) of 
€, = 80z mm-mrad. The electron beam, with a peak current of 140 A is then 
directed into the oscillator (3) where it produces A=10.6um wavelength 
radiation of approximately 100 W average power. [Each electron micropulse 
may have an energy slew of up to 3% across the pulse. Emerging from the 
oscillator, the electron beam is accelerated (4), and then may be compressed in 
two interleaved bunchers (5) to increase peak current. The beam is then 
accelerated further (6) to 34 MeV and the energy slew removed. The electron 
beam enters the power amplifier (7) with a peak current of 450 A and an 
instantaneous energy spread of < 10%. In the amplifier, the optical pulse from 
the oscillator experiences a gain of G=1000 to meet the design criteria of 100 
kW average laser power for 3 minutes [19,20]. 

This chapter examines the goal of passing an electron beam through the 
oscillator (3) without increasing its energy spread beyond 1.0%. This 
procedure produces an optical pulse with 100 W average power, and determines 
an optical pulse shape and size that optimizes the amplifier efficiency. The 
effect of an energy slew in the electron beam is investigated, and we also 


consider the optimum design for a high power oscillator. 


B. OSCILLATOR DESCRIPTION 

The APLE oscillator uses a linearly polarized undulator of wavelength 
A, = 2.36 cm over N=100 periods for a length of L=NiZ,=2.36 m and 
undulator parameter K = eBA,/2amc* =0.23. The optical wavelength 


produced IS determined from the resonance condition as 
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A= Ag(1+K?)/27? ~10.6um. The relatively large optical mode waist, 
W, =0.2cm, compared to the small matched electron beam radius [1], 
= aie. Ay/2nK) =0.08 cm, leads to an averaged filling factor of 
F =(r,/Wo)/(1+ 17/1222) = 0.12 with Rayleigh length Z,=L/2. The electron 
beam density p=14x10'%cm’ gives a dimensionless current of 
j = 8N( neKLJJ)* p/ y?mc* =250 where the Bessel _ function factor 
JJ = Jo(E) — J, (€) = 0.987 with € = K?/2(14+ K?) ~ 0.025 is present in order to 
account for the linearly polarized undulator [1]. The theoretical single-mode gain 
is estimated by G=0.135/F =~ 58%, but is not self consistent and does not 
include the effects of energy slew, beam quality, or optical mode distortion. The 
length of the electron micropulse is /, = 1.8 cm, or o, ~ 17 slippage distances, 
NA. 

The electron phase velocity, v=L|(k +k,)@, —k], measures the resonance 
between the electron, undulator, and optical fields [1]. An electron beam energy 
spread causes a Gaussian spread in z velocities and phase velocities with 
standard deviation Av= 0g = 4aNAy/y =~ 1 where Ay/y= 0.0015 for the APLE 
oscillator. A linear energy slew of 3% over the electron micropulse causes a 
ramp in phase velocities from v= —25 to vx 25. An angular spread of the beam 
causes an exponential spread in the phase velocity characterized by a 1/e width 
of Av, =o, =1.5. The energy and angular spreads do not significantly degrade 
APLE oscillator performance. The large energy slew of Av~50 over o, ~ 17 
slippage distances means that each slippage length of the beam experiences a 


phase velocity shift of Av = 3 during each pass. 


99 


C. DIFFRACTION EFFECTS 

When the effects of optical diffraction are included in the FEL interaction, 
the phase velocity for maximum gain is shifted by the Rayleigh length Zp, [1]. The 
electron phase velocity that gives the optimum peak single pass gain was found 
to be 4" = 2.6+1/z, = 5 [1]. Optical mode distortion results from high current 
which amplifies the radiation overlapping the electron beam. This effect causes 
the optical mode at the end of the undulator to decrease in size and increase the 
filling factor. A single pass, three dimensional, self-consistent numerical 
simulation including diffraction, electron beam radius, and beam quality was 
used to determine the electron phase velocity that produces the peak single- 
pass gain, 4)” , and the new effective filling factor. The simulation used a peak 
current density of /=250, normalized electron beam radius o, = rJnf/La = 0.28 
with a parabolic shape [1], energy spread with standard deviation o, =7 
distributed as a Gaussian, angular spread described by an exponential 
distribution of 1/e width o, =7.5, and a weak optical field. 

The strong fields attained over many passes and the resonator mirrors can 
affect the steady-state distortion that occurs in steady-state. A multiple-pass 
numerical simulation including resonator mirrors and losses was next used to 
determine the optical mode distortion over multiple passes in strong optical 
fields. The losses were characterized by Q=5, where 1/Q represents the 
fractional loss of optical power per pass in the resonator from both outcoupling 
and absorption. The resonator was characterized by a mirror radius so that the 
fundamental mode has edge loss around the mirrors of 1%. Strong optical fields 
cause the gain to decrease until the gain just equals the loss on each pass 
through the undulator. Since the gain is reduced, the amount of mode distortion 


is reduced. The effective filling factor now obtained by comparing the cross 
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sections of the undistorted optical mode with the distorted optical mode obtained 
from the multiple pass simulation gives F = 0.04. It is estimated that transverse 
diffraction and mode distortion in the FEL oscillator can be summarized by using 
an effective jf ~10 in two dimensional codes that explore longitudinal pulse 


effects. 


D. PULSE EFFECTS 

Running the APLE oscillator with a relatively large current of /F =10, and 
extracting only 100W average power, while inducing less than 1% energy 
spread requires the electron-optical coupling to be significantly reduced. The 
first two cases presented below examine options for reducing the coupling. The 
last case increases the coupling to achieve high efficiency without the amplifier. 

1. No Energy Slew in the Undulator 

One method used to reduce electron - optical beam coupling is to operate 

at a large desynchronism, d= —2AS/N/, where AS is the mirror deviation from 
exact synchronism [1]. Small desynchronism values are Known to cause 
unstable operation because a small change in resonator length can alter the 
FEL from high power to low power, or even zero power [1]. Figure 5-2 shows 
optical power versus desynchronism for the APLE oscillator. It was obtained 
using a two dimensional numerical simulation, including slippage effects with an 
electron pulse length of o, ~1/, a Gaussian spread in electron phase 


velocities of o, = 1.5, and outcoupling and losses characterized by Q = 5. 
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Figure 5-2. Power versus desynchronism for the APLE oscillator. 


A larger, more optimum desynchronism of d=0.34 gives low power, a long 
optical pulse, and a small induced energy spread, as shown in Figure 5-3. The 
long parabolic electron pulse of length o, =17 is displayed in the lower-left 
panel showing the slippage over one pass through the undulator from 
dimensionless time r=ct/L=0 to z=1. The middle-left panel shows the optical 
pulse evolution growing from shot noise over n=800 passes to a final peak 
optical amplitude of |a|= 18 shown at the top-left. The shot noise is introduced 
by including a random electron phase with standard deviation of d¢=10°. The 
middle-center panel shows the evolution of the optical power spectrum over 
n= 800 passes with the top-middle panel displaying the final power spectrum. 
The bottom-center panel shows the single-mode gain spectrum for reference. 
The evolution of the electron distribution over the n= 800 passes is shown in the 
middle-right panel with the final distribution at the top-right. The bottom-right 
panel displays the evolution of the total optical power within the oscillator over 
n=800 passes, with the final dimensionless power P=86 corresponding to 


approximately 2200 W. Since the trailing half of the electron pulse never 
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interacts with the optical pulse, most of the energy distribution is kept within a 
narrow band. Only the electrons on the leading edge of the electron pulse give 
up energy to the light which produces the small bump on the left of the final 
electron distribution spectrum. The final energy spread is only about 0.6% and 


satisfies the APLE criteria. 
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Figure 5-3. Optimum desynchronism of d=0.34 produces required power without inducing large 
energy spread. 


The effects of a small random "jitter" in the electron energy coming from 
the accelerator was also explored. The study used the two-dimensional 
simulation but included a small random phase velocity component to each 
electron micropulse. A timing jitter was included by adding a small random 


change to the desynchronism length on each pass. These jitters had a 


=o 


significant effect on small values of desynchronism and caused the FEL to not 
operate for the negative desynchronism values. This effect could be expected 
since the region of small desynchronism in Figure 5-2 is considered unstable. 
The jitters had little or no effect on the power at moderate to large 
desynchronism values. 
2. Energy Slew in the Oscillator 

Another option for the APLE SAMOPA design induces a linear electron 
energy slew of between 1% and 3% along the micropulse. The effect is to 
reduce electron-optical beam coupling in an attempt to leave large portions of 
each micropulse unaffected by interactions with the radiation field. The energy 
slew is converted to a slew in the electron phase velocities by Av= 4zNAy/y. 
The slew in phase velocities over a pulse length of o, = 1/ slippage distances 
results in a shift of 2/4 (1% slew) or 32/4 (3% slew) over each slippage 
distance. Superimposing a Gaussian energy spread of og = 10n the slewed 
pulse causes light to interact with electrons with phase velocity differences of 
~z. Since the change in phase velocity is greater than the gain bandwidth, the 
FEL interaction is affected. In order to reach a stable power and stable optical 
pulse shape when there Is a significant energy slew, the desynchronism must be 
near zero so that the light continues to interact with the same energy component 
of each electron micropulse. Unfortunately, when the desynchronism is small, 
the FEL coupling is large resulting in strong optical fields and a large electron 
energy spread at the end of the interaction. A large desynchronism causes the 
light to experience a significant change in resonance condition as discussed 
above. In that case, a chaotic optical pulse shape and power evolution are 


observed in simulations. 
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Figure 5-4 shows the results of a two-dimensional multimode numerical 
simulation with jF =10, a large desynchronism of d =0.06, energy spread of 
0, =1.9, and an energy slew of 3% over n= 1500 passes. The middle-left panel 
displays the chaotic optical pulse evolution with the final pulse shape shown at 
the top-left. The two upper-center panels show the broad optical power 
spectrum evolution due to a broad range of resonance conditions contained in 
the slewed electron pulse. The lower-right panel displays the chaotic power 
evolution in the resonator with about 50% power variations over a few hundred 


passes. 
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Figure 5-4. APLE oscillator with 3% linear electron energy slew present. 


The three upper-right panels display the electron energy distribution as a phase 
velocity distribution at the end of passes n= 1, 750, and 1500. On the first pass, 


the 3% energy slew and the 0.3% Gaussian energy spread in the electron 
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micropulse are shown. In later passes, the electron-light interaction induces 
two to three times the energy spread originally present, which is the maximum 
allowable by APLE specifications. 

3. High Power Oscillator Experiment 

Another possible goal for the APLE laser system is the operation of the 
oscillator alone at the high average power of 10 kW for 3 minutes. This goal 
requires a complete change in design criteria and objectives from those of the 
previous two sections. Rather reducing the electron-optical beam coupling, the 
goal here is to increase coupling to get maximum power from the FEL oscillator. 

The first method for increasing power simply increases the current density 
j «K* by increasing the undulator parameter K. This is accomplished by 
reducing the gap, g, between the undulator magnets, since K e-9/40 The 
increase in coupling is expressed by / occurs because increasing K will increase 
the transverse wiggling motion of the electrons; the amplitude of transverse 
oscillations is Ax =KA,/2zy. 

Secondly, the FEL power can be increased by using a small value of 
desynchronism as seen in Figure 5-2. 

The two methods above increase FEL coupling and should lead to 
saturation in strong optical fields. Saturation occurs when a significant number 
of electrons have lost enough energy to move through the gain spectrum 
bandwidth Ay/y~-—1/2N, which corresponds to a phase velocity change of 
Av=4nNAy/y~-2zx. Another view of the saturation process uses the FEL 
resonance condition 2 =4,(1 +K?) 23am If a significant number of electrons lose 
enough energy, (y—1)mc*, the optical wavelength no longer remains in 


resonance with the electrons. 
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The resonance condition reveals several ways to go beyond the 
saturation point. As electrons lose energy along the undulator, the undulator's 
wavelength, A,, or the parameter, K, could be decreased so that resonance is 
maintained along the undulator. An undulator designed to take advantage of 
this concept is known as a "tapered undulator’. Conceptually, reducing A, or K 
provides a phase acceleration for the electrons with the effect of simply adding a 
phase velocity to the pendulum equation (2-14) describing the electrons’ 


microscopic motion [1]: 


ie = 6 +|a\cos(¢+ ¢), where ¢ = [ kel eae +kz+at, (5-1) 
0 


where kj(z)=2z/A)(r) is the changing undulator wavenumber along the 


undulator. Such tapering yields a phase acceleration of [1] 


5 eg onal (62) 
az 

In order to increase coupling, the undulator parameter is increased to 
K=0.8. The three-dimensional numerical simulation using K=0.8 yields a new 
effective filling factor resulting in /F = 70. A small value of desynchronism 
gd =0.001 enhances electron-optical beam coupling further. The phase 
acceleration 6=8z continues to extract energy from the electrons after 
saturation. Figure 5-5 displays the results of a two-dimensional numerical 
simulation using losses represented by Q=10 over n=1000 passes. The 
optical pulse is modulated by the trapped-particle instability and grows to a final 
peak amplitude of la| = 570. The trapped-particle instability occurs when 
electrons trapped in phase space oscillate at the synchrotron frequency leading 


to the formation of several sidebands. The final average power P =2.5x 10° 
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corresponds to approximately 20 kW and is twice the high power goal of the 


oscillator. 
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Figure 5-5. High average power from tapered oscillator. 


E. CONCLUSIONS AND FURTHER RESEARCH SUGGESTIONS 
The Boeing/LANL Average Power Laser Experiment (APLE) design required 
limiting the electron-optical coupling in the oscillator in order to reduce the 
energy spread of the electron pulse prior to injection into the amplifier. Several 
methods of limiting the coupling were examined, including the use of large 
desynchronisms and the use of an energy slew in the electron pulse. Large 
desynchronisms were found to meet the requirements of a smooth, relatively 
long optical pulse with less than double the initial energy spread of the electron 


pulse. 
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Using an energy slew to reduce electron-optical coupling resulted in a 
chaotic optical pulse and energy extraction from a large portion of the electron 
pulse. 

Producing a high average power optical pulse using a tapered oscillator 
without the amplifier was examined in this chapter. Further investigation on the 
ability of the mirrors to withstand the high optical powers incident on them, or a 


proposal for a mirror system that can operate at those powers, is needed. 
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VI. WHISTLER:WAVE AMPLIFICATION IN THE 
MAGNETOSPHERE 


A. INTRODUCTION 
1. Description of Whistler Waves 

Whistler waves are very low frequency (VLF) electromagnetic radiation 
that propagates through the earth's ionosphere and magnetosphere, and are 
detected at magnetically conjugate field points on the earth’s surface [21]. 
Whistler waves were first commonly detected during World War | by each side 
as they used amplifiers on long telephone lines to listen to the other's 
conversations [22]. Instead, what they sometimes heard were whistling sounds 
decreasing in frequency from several to one kHz over a period of about one 
second. The long telephone lines were acting as antennas for these whistlers, 
and the high gain devices boosted the otherwise weak signal. It also found that 
the rate of whistler activity varied with time of day, magnetic activity, and location 
on the earth [21]. Whistlers are more common at night, in mid-latitudes, and 
during periods of increased lightning activity [22]. 

The correlation between periods of high lightning activity and whistler 
detection led to the following explanation for their occurrence. A lightning flash 
in the southern hemisphere dumps electromagnetic energy into the atmosphere 
in a short broad-band pulse. This pulse may then propagate through the 
ionosphere and into the magnetosphere. The magnetosphere is_ highly 
dispersive, so as the pulse travels along a magnetic-field-aligned duct it is 
spread out due to higher frequencies traveling faster than lower ones. The wave 


may also be amplified as it propagates along the field line by interactions with 
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energetic charged particles. As the wave emerges from the magnetosphere at 
the conjugate point in the northern hemisphere, it is detected as a whistler 
beginning with higher frequencies and trailing off to lower frequencies. 

Other mechanisms have also been found to generate whistler waves. 
Morse code dashes longer than 150 milliseconds may produce whistlers 
accidentally [22]. VLF transmitters placed in mid-latitudes may stimulate 
whistlers detected at the conjugate point. Over the last forty years, whistler 
waves have primarily been studied using VLF transmitters to simulate the 
injected whistler [22]. 

2. Description of the Magnetosphere 

The immediate environment of the earth consists of a neutral atmosphere 
extending up to about 60 km above the earth's surface, the ionosphere, a region 
consisting of increasingly ionized gas as well as neutral particles extending up to 
about 1000 km in altitude, and the magnetosphere, consisting of ionized 
hydrogen and free electrons, extending out to distances of about 100,000 km 
from the earth [23]. There is no distinct boundary between the ionosphere and 
magnetosphere; it is rather a smooth continuum. 

The structure and processes of the magnetosphere are controlled by the 
earth's magnetic field and the solar winds from the sun. The interaction of the 
solar wind and the earth’s magnetic field forms the bow-shock at approximately 
12R, (where R, = 6370 km is the mean radius of the earth) with the solar wind 
dragging field lines around the earth forming a tail on the anti-solar direction. 

The region in space in which the earth's magnetic field is closely 
approximated by a centered dipole field inclined at 11° to the spin axis is known 
as the inner magnetosphere [23]. The inner magnetosphere extends out to 


about 7R, and is populated by a background ‘cold’ plasma with energies of 0.1 
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to 1 eV and a less dense ‘hot’ plasma with energies from 1 keV to 100 MeV [23]. 
These energetic charged particles are magnetically trapped in helical orbits 
spiraling along field lines and make up the radiation belts. Figure 6-1 shows the 
inner magnetosphere with dipole magnetic field lines and the helical paths of 
energetic trapped particles. The cold plasma is responsible for the index of 
refraction for VLF while the part of the hot plasma that is resonant with the VLF 
wave can amplify the wave. A typical value for the index of refraction is n * 20 


to 30 due to the cold plasma. 





Figure 6-1. Inner magnetosphere with dipole field lines and energetic particle paths. 


Approximating the earth's magnetic field as a centered dipole inclined by 


11° from the axis of rotation yields the following equation for the magnetic field: 
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B(R,A) = a( 2) + sin A), (6-1) 


where B, is the magnetic field at the surface of the earth, R is the geocentric 
distance to the field line, and A is the geomagnetic latitude [23]. A dipole field 
line can be described by 

Rea COSs \. (6-2) 
where F,, is the geocentric distance to the field line at the equator [23]. Field 
lines are uniquely identified by the L parameter 


L= — (6-3) 


The term "ducting" is used to describe the bending of the wave as it 
travels up along a field line in the southern hemisphere and back down the field 
line in the northern hemisphere. Ducts consisting of crests of higher electron 
density, or troughs of lower electron density, nave been found to exist along field 
lines in the magnetosphere [22]. The electron density gradient creates an index 
of refraction gradient within a duct that causes a wave initially entering a duct to 
become trapped much like light in an optical fiber. Since the wave is traveling 
from a higher to a lower index of refraction, a form of total internal reflection may 


take place under the right initial conditions. 


B. WAVE PARTICLE INTERACTIONS 
1. The Cold Background Plasma 
VLF waves propagating through the magnetosphere are affected by the 
hot and cold plasma in the region of interaction. The electrons in the cold 


plasma are several orders of magnitude less energetic than the electrons in the 
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hot plasma and are not resonant with the VLF wave. The only effect of the cold 
plasma on the VLF wave is to cause dispersion. The cold particle density p, is 
several orders of magnitude greater than the density of the hot particles. 

| Since the VLF wave travels along a magnetic field line, a useful 
coordinate system near the geomagnetic equator has the z-axis follow a field 


line, as seen in Figure 6-2. 


PARTICLE 





Figure 6-2. A field line showing the direction of propagation of the wave and electrons and a 
relevant coordinate system. 


The cold particles’ motion is governed by the Lorentz force law 
3S 


mv = a{ En Palen é)) (6-4) 


where ( ; ) =d(..)/dt, q is the charge on the particle, E, =Ex +E,y is the VLF's 
arbitrary electric field, B, is the VLF's arbitrary magnetic field, B =B,2 is the 
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earth's magnetic field, and v is the velocity of the particle. Since the free ions 
are much more massive than the free electrons, they can be considered as 
being immobile at VLF driving frequencies. Assuming the wave oscillates at 


frequency w, the linearized equation (6-4) yields 


~imov, = -e( E, — “By (6-5) 


and 


—IMov, = -E, + “2 | (6-6) 


the VLF wave magnetic field has been dropped since it is always measured to 
be much less than the earth's magnetic field [21]. Solving equations (6-5) and 


(6-6) algebraically for the cold electron velocity, we find 


5 -1 
v= 2 (ie. - 26, 1-3) | (6-7) 
@ O @ 
and 
Ee Qe 1 ay 6-8 
sgt at) Tae. a 


where 2=eB,/mc is the electron gyro-frequency. 
The Lorentz force law shows how the VLF wave affects the motion of the 
cold electrons. The wave equation is required to see how the cold electrons 


affect the wave. Ampere's law and Faraday's law are 








ee eee (6-9) 
Ee Cc ot 
ee 1 6B 
ee eee, 6-10 
ee ca a) 


1a 


and result in the wave equation 


ae eee (6-11) 





-k(k-E,) +kE, = Fi 


C 
where as =-ep.v is the cold electron current density. Substituting the cold 


electron velocity equations (6-7) and (6-8) into the wave equation (6-11) yields 


2 
(aw? ~c?k?)E, =e (5 +i=E sh (6-12) 
1-2 
and 
( 2 22 2) Ope 22 
Os —G k Ey = 7F (6. - SE, ; (6-13) 


where w2, =4z7e%p,/m is the cold electron plasma frequency. Solving 
equations (6-12) and (6-13) simultaneously then yields the dispersion equation 
relating kK and w. From the dispersion equation may be obtained the index of 


refraction caused by the cold electrons, 





Z 
n? = o*k* 2 (6-14) 
w* ola +2) 


The plus sign corresponds to the index of refraction for a left-hand circularly- 
polarized wave; the minus sign corresponds to a right-hand circularly-polarized 
wave. The whistler wave has a frequency less than the electron gyro-frequency, 
and must be aright-hand circularly-polarized wave. 
2. The Whistler Wave Equation 

Having developed the index of refraction, which is interpreted as the 
effect of the cold electrons on the whistler wave, the effect of the hot electrons 
on the whistler wave is now explored. Referring back to Maxwell's equations 


used in the last section, the current density must now include the hot electrons 
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also: J =J, +J,, where J, = —ev,,[-sin@,cos@,0] is the single particle current 
density and © = [alt at is the phase of the electron cyclotron motion. When 
the wave envelope is a narrow band whistler, the slowly varying amplitude and 
phase approximations may be used in Maxwell's equations, QE, << wE,, 
6,E,, << kE,, &¢ << wf, 6,9 << kg, where &, = 2 (..)/6z and G = A (..)/ét [1]. 

Having shown in the last section that whistler waves are right-hand 


circularly-polarized, the fields can be expressed as 


E,, = E,[cos y,-siny,0], By, = B,[sin y,cos y,0], (6-15) 


where y=kz-of+o. Inserting the whistler fields into Ampere's law while 


allowing the amplitude and phase to evolve results in the following equation: 


é,B,,{-cos y,sin y,0|+ B,(k + é¢)[ sin y,cos y,0]- 


1 E 4x - ene) 
— FE wl cos y,sin y,0]+ —(-@ + &¢)[ sin y,cos y,0] = J. 
Defining the two orthogonal unit vectors 
é,=[siny,cosy,0] and é =[cosy,-siny,0| (6-17) 


and projecting equation (6-16) onto them yields the following two scalar, first- 
order differential equations: 
ee 
cB, 2,¢+E, 86+ ckB, - wE, = rj 47ev,,cos(y+@), (6-18) 
= io) 
and 
cé,B, + BE, = 4zev,, sini y+). (6-19) 


These equations can be further simplified by using the index of refraction 


equation (6-14) and nE, =B, to obtain 
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ncé,B, + 6B, = —4zenv,, sin( y+ 0), (6-20) 
and 


B, (ncé,¢+ 6.¢) = —4env,, cos( y+ ©). (6-21) 


We next define a complex whistler field b as 


_ By 
MC 


Then equations (6-20) and (6-21) further simplify to 


b e'* (6-22) 








6 6), _ i4dze? , 
(noZ + Sp - = vexp(-i(¢+ 0)), (6-23) 


where v = nv,, /c is the normalized hot electron transverse velocity and 


C= kz-ot+ {a(t at’ 
is the longitudinal electron phase. 


The single-particle current for the hot electron in equation (6-23) is 
transformed into a beam of electrons by a weighted average over their 
transverse velocities v, their transverse phases 0, and their longitudinal phase 
¢ The averaging is weighted by the hot electron density p,. We also can 
assume the whistler wave is spatially flat over many wavelengths [26] so that 


0 b/Oz ~ 0 to obtain the whistler wave equation 
b = -iw},(v- exp(-i(¢+ ©))). (6-24) 


A whistler wave entering the magnetosphere encounters many electrons 
randomly spread in phase. As the whistler wave equation (6-24) shows, the 
interaction with randomly spread electrons has no effect on the whistler and the 
wave does not grow. However, the whistler wave affects the motion of the hot 


electrons. It is this feedback mechanism that allows electron bunching and wave 
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amplification. The effect of the wave on the motion of the electrons is the next 
topic to be examined. 
3. The Hot Electron Equations of Motion 
The motion of the hot electrons in the presence of the earth's magnetic 


field and the whistler wave fields is governed by the Lorentz force law 


Vp = -£/é, + s x (B,, + é)| (6-25) 


~ However, unlike the cold electron case, the whistler magnetic field is retained to 
describe electron bunching. Since the contributions of the cold electrons has 
been completely contained in the dispersion relation (6-14), the subscripts "“h" 
for hot and ''c" for cold are now dropped after 6-25. 

The microscopic bunching of the electrons due to the wave-particle 
interaction is a small perturbation on the macroscopic gyro-motion and the 
bounce motion of the electrons in the earth's magnetic field. The earth's field 
near the geomagnetic equator is not exactly uniform, because the dipole's 
magnetic field lines diverge as they leave the north pole and converge as they 
enter the south pole. This results in a component of the field in the radial 
direction r, as well as in the previously defined Z direction [24]. In order to find 
the effect of the radial component B, of the earth's magnetic field on the motion 


of the electrons, Maxwell's equation V-B =O is used to find Be 





EB Ce rea == (0 (6-26) 
ror OZ 
which gives 
Ba 4 Se ; (6-27) 
CZ) 5 
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The Lorentz force in the Zz direction due to this field is 


6-28 
z 2 am ( ) 





where v, is the hot electron’s gyro-velocity. Averaging over one gyration and 
identifying r = mv ,c/eB as the Lamor radius of gyration yields [24] 
vi B, 


V,=- —— 6-29 
a (6-29) 


Since the electron energy is roughly constant while spiraling along field lines 


(v," +V a = constant), we have v_v, = —v,v,, so that 


V.V, €B 
Yo es 


= 6-30 
‘3B a (6-30) 
and 
VV 
ened (6-31) 
2B, C2 


Recognizing v~ /(2B,) in equation (6-29) as a form of the first adiabatic invariant 
[24], it is therefore a constant. Thus equation (6-29) can be written in terms of 
the values of the initial transverse velocity v,, and magnetic field, B,, at the 
equaior. 

Adding equations (6-29), (6-30) and (6-31) into (6-25) gives the hot 
electron’s equation of motion including the effects of the whistler wave and the 
earth's inhomogeneous magnetic field. Expressing each velocity component 


yields 


VV 
V, = -Sle, am *(v,B, = v-By)| fee (6-32) 


PBuweze 
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Wz eB, 








e 1 

Vv, =-—|E,+-—-—(v.,B, —v,B,) | + —S 6-33 

5 Sle, +7(v.8, &,)| arr (6-33) 
oh Vio 6B 

v, = -—-\v,B, -v,B, J -— —. 6-34 

enc mee” 2b oz aun 


Since the region of space most often identified with the whistler wave- 
particle interaction is within a few degrees of the geomagnetic equator, a useful 


approximation to the earth's dipole magnetic field in that region is 
B, = B(z) = B)(1-qz?), (6-35) 


where g = 9/ (21?R,?), B, is the magnitude of the earth's magnetic field on the 
field line at the equator, and z is distance measured along a field line from the 
equator [25]. 

Substituting the form in equation (6-35), for B, as well as for B, and B, 
and the right-handed circularly-polarized whistler wave from equation (6-15) into 


equations (6-32), (6-33), and (6-34) yields 


Ne ie -£{ Ecos y+ “Y B(z) ~ “2 nEcos v| +V,V, QZ, (6-36) 
m C C : 

S . V V . 

= -—| -—Esinw-—B(z)-—+nEsin VV, QZ, 6-37 
Vy al =e (z) ; vs y¥z,9 (6-37) 
: e;v V 
Vv. = -—| —nEcosy—-—nEsin y|-|v7_}qz. 6-38 
2 <(t ae: v|-(i)a (6-38) 


These equations can be further simplified in terms of the perpendicular velocity 


v, and the angle © between the whistler wave fields and the electron’s velocity, 
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vy, =v, cosO, Vv, =v, sino. (6-39) 











where © = a+ {al ‘\dt'. After some algebra, the equations of motion reduce 
to 
° o. Carva 5 
V, = —-Vi qz+—*“+ sin y+ 0), (6-40) 
: mc 
eB nv, 
= “| 1- siniyw+0 6-41 
avig2 + Bu (4- Me sin y +0), (6-41) 
- @eBlz 
6 = SAE), 2Be (4 Mz M2 Joos( y+ 6). (6-42) 
me nmv, c 


Introducing the whistler wave amplitude |b}=eB,/mc and the electron gyro- 


frequency Q = eB(z)/mc yields 


V, =-Vvi gz +|blv, sin(y+®), (6-43) 
Nn 

VE = Vevinge AI s- — }si y+ ©), (6-44) 

6 = A(z) +L! 1" S loos vee (6-45) 
nv, c 


The first terms on the right side of each equation (6-43) through (6-45) describe 
the fast macroscopic electron motion in the large scale dipole field of the earth. 
The last terms are proportional to the VLF wave field \b| and describe the small 
microscopic scale motion in the combined dipole and VLF fields. 

Since the hot electrons spend much less time in the interaction region 
than the bounce period in the dipole field, the electron motion in the first term on 
the right side of equations (6-43) through (6-45) can be described by 

Zmacro(t) av, t. (6-46) 
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According to experimental observations, most whistler activity and VLF 
wave amplification occurs when w =0.52, known as the "cutoff frequency” [22]. 
At the cutoff frequency, the phase and group velocities of the wave are given by 
V,, =¢/n and v, =dw/dk=c/n respectively. Since at resonance v=0O, the 
electrons satisfy kv, = @—Q so that v, =-a@/k=-v, =-c/n. Hence, the wave 
at the cutoff frequency and the resonant electrons are traveling with equal 
velocities in opposite directions. 

Recognizing that © has both fast and a slowly rotating terms allows us to 
separate out the simple, fast periodic gyro-motion of the electrons in order to 
concentrate on the more interesting slower evolution of the bunching electrons. 
t 


Redefining © = 6, +6,, with 6 = | 


Olt )at’, allows the use of 6,(t) for the slow 


evolution of the VLF wave. 
The arguments within the trigonometric functions in equations (6-43) 


through (6-45) can be rewritten in terms of the longitudinal electron phase 
C= kz~ at + f O(t")at’, (6-47) 


where y+@O = y+6+6,=¢€+6,+ ¢. Defining the electron longitudinal phase 


velocity as 





; 242 
=e 0+ 04 14 x | (6-48) 


where (6-46) is used to replace z in terms of tf in the inhomogeneous term. Thus, 


the electron phase acceleration becomes 


oe Oo 
C= v= kv, + 4q7t. (6-49) 


no 


Substituting the electron phase and electron phase velocity into equations 
(6-43) through (6-45) while using (6-46) for the macroscopic z terms, v,, = ao|k 


for the macroscopic terms, and v, = -—c/n in the microscopic terms yields 





v= (4- vi) +|blva sin(f+ 6+ ¢), (6-50) 
2 

y= NO t+ Qlblsin(l+0+¢), (6-51) 

6= Al cos(¢+ 6+ ¢), (6-52) 


where the definition v= kv, is a normalized perpendicular velocity, vo = kv is 
the normalized perpendicular velocity at the equator, and the subscript s is 
dropped from 6,, since all coordinates are now slowly evolving. The 
macroscopic z motion along the inhomogeneous field line uses the substitution 
z=ct/n. The microscopic z describing electron bunching in the VLF wave uses 
the fact that the wave and electrons are counter-propagating which causes 
Vig = Ole 

The first term on the right side of equations (6-50) and (6-51) represents 
the acceleration of the hot electrons due to the inhomogeneity of the earth's 
magnetic field. The acceleration term is analogous to the “taper' term in a free 
electron laser as discussed in Chapter V, and so the acceleration term is 


referred to as a"'taper" here. Define 


2 Z 
aq Io 
=- Aan 6-53 
kono sians oe 


where L defines the number of geocentric earth radii at the equator to the 


magnetic field line being used, and R, =6370 km is the mean radius of the 
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earth. Presenting the hot electron equations of motion using the taper term 


gives 
v= (4- v2) dat +|blvasin(¢+ 0+ ¢) (6-54) 
V = vot + 2\b|sin(¢+ @+ ¢) (6-55) 
6 = Al cos(¢+ 6+ ¢). (6-56) 


Several interesting observations can be made about the equations of 
motion (6-54) through (6-56). Equation (6-54) follows the longitudinal evolution 
of the electrons, and includes the pendulum equation with a phase acceleration 
identified as taper. The taper term in equation (6-54) shows that certain 
electrons with an equatorial normalized perpendicular velocity of v)g+2 
experience no effect on their longitudinal phase velocity due to the 
inhomogeneity of the earth's magnetic field. A normalized perpendicular velocity 
of Vp = +2 corresponds to a pitch angle of = +63°. The electron phase velocity, 
v= kv,(t)-— @+Q(t), shows that the parallel velocity and the gyro-frequency 
are each changing in opposite directions at the same rate at this particular pitch 
angle, thereby canceling. The perpendicular velocity is only slightly affected by 
taper since (6-55) evolves much slower than the longitudinal phase velocity 


equation (6-46) due to the presence of the large factor w = 10°/s. 


C. NUMERICAL SIMULATIONS 
The equations of motion (6-54) through (6-56), together with the wave 


equation (6-24), are now numerically integrated. The solutions reveal many of 
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the same characteristics found in free electron lasers, including the observation 
of electron-bunching, saturation, trapped particle instability, and taper. 

The parameters used for the simulations are taken from measurements 
made by the Space, Telecommunications, and Radioscience Laboratory at 
Stanford University from Siple Station, Antarctica. This data was obtained using 
their >100 KW transmitter to send VLF signals at between one and six kilohertz 
using a long horizontal dipole antenna [26]. The amplified signals were then 
received at the conjugate point at Lake Mistissini, Quebec [26]. The data being 
simulated here is a transmitted signal at w = 10*/s estimated by the observers 
through time of flight measurements as traveling along approximately the L=5 
field line. An index of refraction within the interaction region close to the equator 
was estimated to be approximately n=25. The hot electrons present are 
represented by electrons uniformly spread in the electron phase, ¢, =0—- 2z, 
and the electron cyclotron phase, 6, = O-—> 2z. The time of interaction is taken 
to be 7=1s. 

The first case examined is in the weak-field, high-gain regime. In an FEL, a 
weak optical field is identified by lal <7z. Similarly, a weak whistler field is 
identified when |b|w7* < z which corresponds to |b| ~ 3 x 10~*/s [27]. Figure 6-3 
shows the results of high gain in a weak whistler field of initial value 
by = |b(0)|=5x10-°/s as it interacts with a hot plasma of frequency of 
@, = 0.2/s. The electron beam has a small Gaussian spread in phase velocities 
characterized by width standard deviation o, =0.5 /s about the initial phase 
velocity vy, =O. All electrons have a normalized perpendicular velocity of v, = 1 
corresponding to an initial pitch angle of 2/4. The whistler’s frequency is 
@ = 10° /s with an interaction time of T=1 second. The taper term is 5=0, 


which indicates a homogeneous magnetic field. An FEL, with current density 
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j2 7 is in the high gain regime and experiences exponential gain. Similarly, 
when @7aT° > z, (oF a, > 7/@T° = 0.02/s) the whistler wave is said to be in 
the high gain regime where exponential growth occurs [27]. After a characteristic 
time elapses allowing the electrons to bunch (approximately 0.35s in this 
example) the gain becomes exponential. In an FEL the exponential gain is 
calculated from = ¢¢, (z) ~ exp ( j/2)"" V3z|/9 [1]. Similarly, the exponential 


gain for the whistler wave is calculated by [27] 


G(T) = exp| (wo /2)* ar / 9. (6-57) 


Equation (6-57) yields a final gain after the one second interaction of G = 34 dB. 
This value is slightly larger than that shown in Figure 6-3 since equation (6-57) 
does not take into account gain degradation due to a spread in electron 
longitudinal phase velocities or the delay due to the bunching time. The phase 
is shifted linearly by nearly the same factor used in the exponential to calculate 
the gain, Ag = (ow? ID ey 2 ~ 3 [27]. The phase space in (¢, v) displays the 
final positions of 2000 sample electrons after the one-second interaction. Those 
electrons going up in phase velocity remove energy from the wave, while the 
electrons going down in phase velocity transfer more energy to the wave and 
Cause amplification. Tne bunching of electrons in phase is clearly present at two 
locations in phase space since there are two phase factors (¢, =O—2z and 
6, =0— 2z) within the trigonometric functions in the equations of motion. The 
two dark horizontal lines on the phase-space plot show the maximum height of 


the closed-orbit region given by Av, = 4,|\blaT? (see Figure 6-3). 
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10 log (1+G) 


Figure 6-3. Whistler interaction in a weak VLF field with high gain. 


Most whistlers reach strong field saturation during the interaction time. 
Once saturation occurs the electrons become trapped in closed orbits in phase 
space and execute synchrotron oscillations at a frequency , ~ Job] [27]. 
These oscillations can evidence themselves as oscillations in the power and 
gain of the whistler wave. Figure 6-4 shows the results of a Gaussian 
distribution of hot electrons of width o, =10 /s about the resonant phase 
velocity v, =O interacting with a strong whistler field for a time of T=0.4 s. The 
electrons are at a pitch angle of 2/4 determined by v,=1.0, with a density 
corresponding to a hot plasma frequency of w, =10/s. The whistler wave has 
an initial amplitude b, = |b(0)|=0.01/s with frequency @ = 10° /s . Again, there is 
no taper so 6=0. The phase space plot of Figure 6-4 shows the final positions 
of 3000 sample electrons showing their bunching at two phase regions and each 


with lower energy. The spread in phase velocities causes the bunching to be 
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more diffuse and creates much of the random scatter of electrons seen in the 
phase-space plot. Electrons appearing outside the final separatrix lines have 
been trapped by the previously larger whistler field in the oscillating pattern seen 
in the gain plot. Note that the final gain is at a minimum caused by the 


synchrotron oscillations. 


xx*x*x Whistler Phase Space **** 
bo=0.01/s W=10000/s ven Vo=l 
T=0.4s @,=10/s 5=0 0, =10/s 


10 log(1+G) 





Figure 6-4. Whistler interaction in a strong VLF field with high gain. 


The beginning of the gain curve shows that after the short bunching time, the 
gain is exponential until strong field saturation occurs at a value of 
Ib], = 2( wi ao = 1.3 /s. The strong fields trap the electrons and cause them 
to perform synchrotron oscillations. The presence of approximately 4.5 
synchrotron oscillations is clear in the gain plot and agrees closely with the 


calculated value of w, ~ ./@|b| ~ 36 /s which would predict 7 full oscillations 
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over the interaction time of 0.4s. The small discrepancy between the number of 
synchrotron oscillations observed and the estimate @, ~ 36z /s occurs because 
of the spread in electron phase velocities. Saturation in strong fields, electron 
trapping and synchrotron oscillations are all nonlinear effects observed in both 
FELs and whistler waves [27]. 

The previous figures used a short time of interaction and a homogeneous 
magnetic field with 6=0. For a longer interaction time, the effects of the 
inhomogeneity in the earth's magnetic field acts as a phase acceleration 
analogous to taper in an FEL. In Figure 6-5, the electrons begin away from the 
geomagnetic equator and far off resonance with a Gaussian spread of phase 
velocities characterized by a width o, =10 /s centered at v,=2360 /s. The 
electrons have an equatorial pitch angle of 2/4 determined by v, =1, with a 
density creating a hot plasma frequency of w, =10/s. These hot electrons 
interacting with an initially weak whistler field of b, = b(0)| = 0.01 /s of frequency 
@=10* Is for atime of T= 1s. The taper value of 6=0.27 corresponding to 
the duct at L=5 with index of refraction n=25. The phase-space diagram 
displays the final positions of 3000 sample electrons. The electrons begin far 
above resonance, but their phase velocity decreases as they move into 
resonance near the geomagnetic equator. The simulation is designed so the 
sampled electrons are in resonance with the wave at the geomagnetic equator at 
atime t=7/2. After strong coupling between the hot electrons and the VLF 
wave near the equator, the taper causes the electrons to increase their phase 
velocity as they move past the equator and away from resonance. However, as 
seen in Figure 6-5 most of the electrons are trapped by the strong whistler wave 


field and continue to give energy to the wave. 
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Figure 6-5. Whistler interaction in a inhomogeneous magnetic field with a strong VLF field and 
high gain. 


The gain evolution in this figure shows no coupling between the electrons and 
whistler until the electrons’ phase velocity is nearly resonant at time f=0.5 s. 
Once sufficiently near resonance, the electrons become bunched and the gain 
becomes exponential. Synchrotron oscillations begin. The effect of the taper 
causes the trapped electrons to continue to give up energy to the wave after the 
normal saturation. Most of the electrons remain trapped near resonance by the 
strong whistler field; even though the earth's dipole field provides a phase 
acceleration to move them away from resonance. Those electrons not trapped 
by the strong VLF wave are moved off resonance by the taper. The gain 
evolution shows that because the electrons are far off resonance initially, there 
is no growth in the wave. As the electrons approach resonance at the equator, 


the same exponential growth seen in the homogeneous case occurs in the wave, 
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leading to saturation and synchrotron oscillations. The inhomogeneity of the 
earth's magnetic field provides a phase acceleration that allows continued 
energy extraction from the trapped electrons. This is seen as the gain slowly 
increases even after the homogeneous saturation level. The wave's phase 
initially decreases, a characteristic of tapering, and then begins to grow quickly. 
The rate of phase growth decreases in conjunction with decreasing gain growth. 
The general characteristics of the gain and phase evolutions shown in Figure 6- 
5 are in agreement with those obtained by Carlson et. al. of Stanford University 
[28]. Their simulations also show gain initially growing exponentially and then 
reaching either a steady value or growing slowly [28]. They also found the phase 
would initially decrease during the exponential growth, and then would increase 


parabolically [28]. 


D. CONCLUSIONS AND FURTHER RESEARCH SUGGESTIONS 
Although apparently an unrelated process, the amplification of whistler 
waves Is shown in this chapter to be closely related to the amplification of an 
Optical wave in a free electron laser. The radiation field in both whistlers and 
FELs couples in resonance to "wiggling" electrons for efficient energy transfer. 
Nonlinear effects such as saturation and the trapped-particle instability are 
observed in both. For the whistler, the inhomogeneity of the earth's magnetic 
field lines acts in an analogous manner to tapering an undulator. The electron 
equations of motion and the wave equation derived in the chapter are used show 
many of the same effects seen in whistler and FEL research. The results closely 
reproduce the results of other researchers using more traditional approaches to 


the study of whistler wave amplification. 
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This chapter solves for whistler wave amplification of a single mode using 
either monoenergetic hot electrons or electrons in a Gaussian distribution of 
phase velocities with a single pitch angle. Further, all electrons are forced to be 
in resonance at the equator, or at time tf = 7/2, regardless of their initial phase 
velocity. The research needs to be extended by adding a more realistic 
distribution of hot electron energies and pitch angles. A multimode numerical 
simulation should be developed, similar to those used in FEL research, so that 
the experimentally observed sweeps in frequency of the whistler waves can be 


investigated numerically. 


89 


[1] 
[2] 
[3] 
[4] 
[9] 


[5] 


[7] 
[8] 


[9] 


LIST OF REFERENCES 


W.B. Colson, Chapt.5 in: Free Electron Laser Handbook, ed. W.B. Colson, 
C. Pellegrini, and A. Renieri, North Holland Publishers, 1989. 

J.M.J. Madey, J. Appl. Phys. 42, 1971, 1906. 

H. Motz, J. Appl. Phys. 22,1951, 527. 

R.M. Phillips, IRE Trans. Electron Devices 7,1960, 231. 

L.R. Elias, W.M. Fairbanks, J.M.J. Madey, H.A. Schwettman, T.I. Smith, 
Phys. Rev. Lett. 36, 1976, 717. 

D.A.G. Deacon, L.R. Elias, J.M.J. Madey, G.J. Ramian, H.A. Schwettman, 
and T.|. Smith, Phys. Rev. Lett. 38, 1977, 892. 

W.B. Colson, Phys. Lett. 64A, 1977, 190. 

W. B. Colson and J. L. Richardson, "Multimode Theory of Free Electron 
Laser Oscillators", Phys. Rev. Lett. 50 1983, 1050. 

W.B. Colson, in: Free Electron Lasers: Critical Review of Technology, ed. 


B.E. Newnam, Proc. SPIE 738 1988. 


[10] J.D. Jackson, chapt.14 in: Classical Electrodynamics, John Wiley 


Publishers, 2d ed. 1975. 


[11] W.B. Colson, Nucl. Instr. & Methods in Phys. Res. A296, 1990 348-350. 


[12] G. Dattoli, Chapt.7 in: "Free Electron Laser Handbook", ed. W.B. Colson, C. 


Pellegrini, and A. Renieri, North Holland Publishers, 1989. 


[13] P. Elleaume, IEEE J. Quantum Electron. QE-21, 1985, 1012-1022. 
[14] G. Dattoli, A. Renieri, A. Torre, J. Gallardo, Phys. Rev. A35, 1987, 4175- 


4180. 


[15] G. Dattoli, T. Hermsen, L. Mezi, A.Torre, J.C. Gallardo, Nucl. Instrum. 


Methods, A272, 1988, 351-363. 


90 


[16] G. Dattoli, T. Hermsen, A. Renieri , A.Torre, J.C. Gallardo, Phys. Rev. A37, 
1988, 4326-4333. 

[17] G. Dattoli, T. Hermsen, A. Renieri , L. Mezi, A.Torre, J.C. Gallardo, Phys. 
Rev. A37, 1988, 4134-4139. 

[18] A Zauderer, Partial Differential Equations of Applied Mathematics, 2nd ed., 
John Wiley and Sons, New York, 1989, pp 48-54. 

[19] D.C. Quimby, C.G. Parazzoli, J.B. Romero, B.D. McVey, and D.A. Byrd, 
Nucl. Inst. and Meth. A318 1992 680-690. 

[20] D.C. Quimby, C.G. Parazzoli, and D.J. Pistoresi, Nucl. Inst. and Meth. A318 
1992, 628-635. 

[21] L.R.O. Storey, " An investigation of Whistling Atmospheres", Philosophical 
Transactions of the Royal Society of London, 246, 1953, 113-141. 

[22] R.A. Helliwell, Whistlers and Related Phenomena, Stanford University 
Press, 1965. 

[23] U. Inan, “Non-linear Gyroresonant Interaction of Energetic Particles and 
Coherent VLF Waves in the Magnetosphere", PhD Dissertation, Stanford 
University, California, 1977. 

[24] F.F. Chen, Plasma Physics and Controlled Fusion, Volume 1: Plasma 
Physics, 2d ed., Plenum Press, 1990. 

[25] U.S. Inan, R.A. Helliwell, W.S. Kurth, “Terrestrial Versus Jovian VLF 
Chorus; A Comparative Study", Journ. of Geophys. Res., 88, 1983, 6171- 
6180. 

[26] R.A. Helliwell, "VLF Wave Stimulation Experiments in the Magnetosphere 
from Siple Station, Antarctica", Reviews of Geophysics, 26, no.3, 1988, 551- 
Se. 


ail 


[27] W. B. Colson, B. Gately, D. L. Caudle, and K. A. Sturgess, "Analogies 
Between the Free Electron Laser and Whistler Mode Amplification in the 
Magnetosphere", presented at the 14th International Free Electron Laser 
Conference in Kobe, Japan, August 21, 1992. 

[28] C.R. Carlson, R.A. Helliwell, and U.S. Inan, "Space-Time Evolution of 
Whistler Mode Wave Growth in the Magnetosphere", 1990, 15073-15089. 


92 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 
Cameron Station 
Alexandria, Virginia 22304-6145 


Library, Code 52 
Naval Postgraduate School 
Monterey, California 93943-5002 


Professor William B. Colson, Code PH/Cw 
Department of Physics 

Naval Postgraduate School 

Monterey, California 93943-5000 


Professor K. E. Woehler, Code PH/Wh 
Chairman, Department of Physics 
Naval Postgraduate School 

Monterey, California 93943-5000 


Professor John Neighbours, Code PH/Nb 
Department of Physics 

Naval Postgraduate School 

Monterey, California 93943-5000 


Professor Robert Armstead, Code PH/Ar 
Department of Physics 

Naval Postgraduate School 

Monterey, California 93943-5000 


Professor Rudulf Panholzer, Code SS/Pz 
Department of Space Science 

Naval Postgraduate School 

Monterey, California 93943-5000 


93 


8. Professor Maurice Weir, Code MA/We 
Department of Mathematics 
Naval Postgraduate School 
Monterey, California 93943-5000 


94 








DUDLEY KNOX LIBRARY 
NAVAL PO 7G ADUATE SCHOOL 
MONTEREY CA 93943-5101 



























SUM yen) 
J at | Sa 
af Pt be 
ef Bh bal 








8. ge 
te eee 

re Ps 
a ae es pon 



















eed | 


xa 3) Siar, a] 


i aM cre at 4 
pact seed hy 
























' NNN clk | 


3 2768 00307424 6 













FI 
(iC ia 
~ ‘ a ah b> « 
ad ROR Tuh a hee a) eee Neem aan ot try mi 
Deas ANS het ? - 
tabate ge ae bk ' » 
Sento . 


) 
Pepe liseiah 
ot rw cepa 
iets ee yet 


cas 


























Saar ie rf wy 
4 1 MD) oe Loe Peed ee | 
oe ie rf ie cers 






es 


ale Oe areas 
a 
ey 


its tren ee, Pony, to Ry, 







































t 
U Lacey f 1 et 2 Fi ed ek 
ne PE be eT CY ER, EPL re ey 

MRD a cel Uae ole tielne * 


ris 






more ee Tr 
aT el 
Tat F xy 


rats taeae 














‘epee fee 
dr’ 4) oe 
CET aTCy "4 ae any ae 
1 Lt , 
Hd A ar a PRET Red 
5 she. 5 7 
NAN Wave PAT ewer Lat Se eas 1 Pop e 










Oe oe 





nT Pi i 
ee RAF oa biti iberrnaees rn 





































































































a 
A ¥ ay rr Le Pn 
7 rer tte oe te a4 
ae ators PAV coil CREMP 
W a ne 7 2 bd 2 8 Gt a ¢hee we + L * 
The Are / Cre j D eas ear ee Ucpaad Wy eA rs ne, 1 Coe 
fe PAs ved Pr ry oy Ls OA ti tbh aa 4 ry ener) F 
oP a See rn y pte yA " an) A : 
« ay) ee Red Ae Sree AY So Me Pa ls ie 
Ltd HE PH 7 ae eit pettea 
DsFee Pareto] yf hae EP Ae 
obs 3 198 Sad De fe E P A 
Af 1 bd is AJ 
tee mike ty Y iehs UF Mey " if wae Rath + cH ty 
Sept Seeley i Artes ries ee Ca 
A Y iy ‘ ( 
oo te de ay as rae RIA Trt pag ee 
I 7 on ’ 0 ’ = cl 
Pier oe oe ae tae 4 re re ea Lp 


er aa 
Pha + i} 
ear ie ‘; Pi : re a 


$i He 












WT aaa ye erry ty ib H 
tt: fhe ie; PUD Se) oe ee a £7, mrt D ta 





















































































































































































































































eat dee ae oe 
Past ae 


Cee SF ne 


ga 
ae ca ry rh Pe nf a 7 Ad he 
hy ee ut Seige ett Es oe 7 ro 
oO vo a ft +? roe ry FP 
ie a 1H ez ASD HONE Sete sent ran , oI ri Ps oY Fy oe coe ‘ th a 
mits A Peiyet NGA BER at rT iy mY ae Bf Cit asees yx bd testtassatate$ ru ear ere sdebals 
pepe HAS oe aM ¥ ue = Ba ‘ ry erate of bets gud eee urn vi PA ea ORLY allie heat eA aie 
poe ree Lean ; y shed af io ab palo ue er ALTE seme Barer oa ee are ara 
Srey ney Y oo ee pay Veh, neat ds dak aif H Saath or an: ane T Tara tO Oe ed ere ora. Ee et ek ee 
t Ppt are ah is lite aga: Fe ese eM Hadee et AE es tos yn | 
eed La daregy Tere % _ Poe ee Oe Ie a ae ct ee 
» Lae Aid es on ine ay mo #7] ne RUS Le J bee , 
4 f ! a ieer tt aris onde ree oom 
Be etait oA. Sat a at os eit Rbk vk ee Oe yh ie, ahr perce peat 
; ae ie Aue b t hg at Wa Rea | a mre AR Sree 
fh ae ne ws a mkeh Ry Ey it elk Lae gee eh POW FFE abe ra Roark rents | Cg Ha 4 
7) th ae Pa Pre ty aS Ra a ae ee sa eS Us orent anh Mens , te Rae) cay i 
ae rf a Ate tere ‘ress hee i ai mn : eS OS aN Bt f Fo cae a ts col Pores 
x f 2% . Cy, 2 ic A Ps 
Ly ide td aaa eagantt sagged st se Bes pirttd os i eRe any “Cra ie of fa ae eee Ag ae te a 
Dati Mbt an mt ‘i BON ees aL oi Tiree iu Re an Ay , Sateen rma 
- ry o braleiata isin, et P +) i Pom biege te rH aos fries Hance a ne ie ere ‘ bs 
Pah hoes i? PPt tel Jud Fe asi is Hee sant A " ar Arpt eae ot 
set th x. Cty a hey Fisch. oe bit tre A a ae if a | a Sr a.) Ly eae CR 
Pe ve ae ‘LP ed gam Pee ar ac ne lepeteye ae 
abbas nae . i Pe Chepy ed Od © or Oe ; Ltn’ ee pe a i : 
re DV ele pedr aru opy Lae a rt i 3 
i aia as aca Nea to ; 
es aie Lia eT} fe ’ oy od se i 


ay cme ee Oe i} 


oo tgs M 
ees yh paeas et a is , 











bit 
“begs plesa® 
ve a tg) 











tal 
ase iii ete iD 
eet PRPS Peete) ha 
Le iele tht ee mp 
Pa Eey i IaerPtiur pn! 

aa ae 28 Ceres ae Sees 


} a 
PEE ALE eee 



















As 
od ae 
fei ncaPetar'y ere ty ie 
a4 nae, ke Se Ye 
ee er ee 


A (eae 
PGS a | 



























aie boy HAN 
































































































rere tate er eh wey a ee 
sO Fr Gut he sti ran Sates rw Te on hye Tre meed ce i 
3 & ie ery "Aky Aa ae ‘ob fee ener yy i ; part 
Ores ay Hebd Uy ar TS YTe rhe Hated Meh Ra en Tt ae a ts te CRD To old ' A Ue Bees ore Ab A 
PPA itt oR Theda hatceee Pee ary oe are a pe HY Nt Nate yNeba de 
Ue Pe ee ee Ae TT) tg Tee Tae) 5 rT os CY Ld j Ps eA an Re CL atom 
es ttre) bast tone by ib ea Rane Oc 
SOS heb nates Me bd Fa ee Te uJ E AM sees ase ot J 
Rh tarts yeast 14) Lae en Pe eae no a a 
cn eee pist rck bor Oe aad yh ae arias a erat } 
oe FEY Lert? ya sats oe ry Hy irae 
4 3s rk (ot tia Yd ne ee Ri 
eo Pig tats § ee eats - 
FT 4 bade he ae eat ahh fe Fd ee re rth i 
* ns 
“a cg 
- Se 
Tad 


PLO nef age ort ee 
Pee date ETO: 


Ch rea ty iat, 
r, NEA Bt See 





























rn ‘Yea rh Ee? ee to 
a) hela be ee ane Pe 4 Raa wy 
H ad ee ran ‘ai 1 ve vy Were Shree a hoe as H 
A ye 2 Carre ey Seu re Fes 
ah i Wn aL Wires hy iY wide Ar Be Lae 
Pot At ae mT 


ETE OM Ad a A att ae He 
eri Sa os Yat make. 7 Pyrarsy 


a TENT Soe) Lae 
é hate ae by He 


ere a a 
data er ae eae eS te Ce Ree 
par A ae Fes b 

re ; 









han 
Cd ee Seas bee 
Ys apt ee ga rate 
he ees td ae “hy 


oes 
















“ 
1 
Le nS fae 
be 






















































































































a 
ee acaba, HE stee foe rae herr patil 
F a et aie ER ae? OT ae Mi ry hb nthe OT aie 
ra og! i Yea era Py a \ A A or 
KE Pana yr La oe hes Ere eee bee, A! us “ } A far rier! al) A sa at ap 
Sr ae eater ro rey by : ; ahd Ble is a Vega gta hae Pest ys dae, 
itera ES nak ney hf UAE HA Meera eros ea a A PG or bt 4 So as dlp ae Nore i an) oer 

cf ; Hd da be LT Oe YN Bs p eth f LR a gil CN 
ayia AS eee 4 ee eo pare he fy AYES) a a Miata aie 4 es Sate HP aera pe gst Rae ea 
Misr . ht at », aay Ca aie iy ae 7 rs BPRS Ane Reeve Onde) “> be’ teeh babae 
arta x bay oy H i 


rh PLoS Sees Ah 


i ar oo 
3 Pn hat) at « 
ry 


» ae a 
eA Te Ret a a eee 





a ee at 
ees ta Ct 





oer ee ee 
Sri 


eee LT oy 





















Pt eee et 






















































































































































































































































































































































































































































































































































ay dh a 












eerie Sa fies « rh 
F : a peed weg Ms ° > See 
: ar Ce oe ' 
ra PS es Say ie ee ory 4 e rs 
“75 eeie at Ey PeNs yeas toe, Bitty 5 phe aa ee Ag 2 
WP MIN oo, oe nee 4 Ste 2 ee gs b , fa = He ts 
Ah ed 9 roars t a a ote & Pie - a BA on 5 elite s 7 
re ed be | » beet od, A ae Me magnet Ae Saye U A ates 8 5 ed] ’ 
ur SACLE UY, ot rae fa 8s 2 Tae ee ee ve < Prim it? ai & pr 
Bp ede rey s pre eerie oars Pt) Cm HL far : : : is 
y ere i 7 rae tir ‘ yes om teres Lb Fi 
te EA Tee Pee 1 Wirt Foon adamant d ¢ S : 
beans aad = At . rT 7 z f 8 a 
‘ WP TNs Meet oe anh Phat Sdees Pres 
ha ba Paes et Bo \e beg ange tere ries re Fee Pe aa} =a ‘ Y 
a zoo ot Coes Th Pen IS he * yr ays arr 7 weet oan PP 
stay AE Ter Ai rol .? ra G ; 
, f ron > Pi ’ by 
1 Wicd Pa tera Vit Har SoBe AN mee tatpey  agt ir 
nT Pier, ee reeverp g a H ailpe Pies ae * 
Pe, peas Par) Bane er) aS eae SP weet) i ee a er rth eres Mart | oe 
ry 4 . ce 1s , ‘ m4 Pr ’ “ 
er ak ae ion a cdl es EE RAE a pa : i Ne Pha Hel 4 ey ee Rar Hi = s * . = ; 
; ; A as 
ya eres ou eek P et SUR Ser nan aR OE ca Br Litt. rary) (a0) tens mk , 7 mr 
Sy 4 of, 7) erties ad a eS an ee 1° wrt ae ‘ the b 7 * ry 
a atte Hy he oe Le ry ard rf ; ee Fd iv ' 
ry et oh ey S rar “spacxbr ey ek lb 4) ti el he ee o i 3 
MH UO stata eee & tg ae One his r 
tt Lay ie bs fe eT earin t ar Ae et v7 ty cree 4 Pore ND ‘ Pa 
ae 4 es Bs a3 yb mute Rete xs, ER A! tesa ls set? so ge fay 8 Te ’ ay . 
Fs ' "ts " P ; F G - 
to oi hd FN) ‘ a es #, othe es neh EL ea Se erp oe ; peo n or in 4 i “ anh nee A roan Ca 
ae eee SP Ts tT eee T re! bob hh ce eee tte SAT ae yd tre fd ry © & pent i 1 " Teen i re re. rt: HM y 
i eta ded BAe ee ead ae beh oh Aby ed ea wy a" tf os = . , Pr tae Py: ha] ? 
OLE Te ae cree a ots LTE Eh ce en oe ed cn Gan eens ee APA try ere Le = 
ed Sea cA A a ee ee as arotlan a oe EES ry ra ee é a fm eRe Pv 
areyri tap ay a hare a) 4p} Ce] ae) , ’ ‘ FY Uae Tac = 
ri ‘ ° r te rs af ’ rt ‘ be s s 
} | fet Oe dee | a) Cn | 7a 
my > 7 ar ie ad 7 My 
ee eat atite Py es i et Fy a] 
eae : an 
A air etree fi rT Ta o ee ah | 
eth at Rarer wt Cw ae atts Py toner ar P 7. a 
Sy Hyer LOEW eter er H — Fi } . 
ee at WPM LPR try eee i Rey ing as) “f ame Ca PO a ea HE , P i é r . a? oe 
neh Bo ae ear ee ape mh eed es , EPL ae Se oe a a Be aa cae i eCine ey oye “un + Ca ae Co a 
a ‘ F 5 - 
Suse es BD oF a Seery) ar a Te EM ean eH sea 7 Un eea To Ay RUS “ 
i ate eat Reed i thet ed as oe pees he Pete ok iy beet aae my ia * ® Oo 
pag a7) Les Sh! * rare Ps iat . , 2 
AN F398 8: 8 ae AS a He Pe Ce ect cS ait sth. or + : Fd ete ar Meo 
a eater ra ae Weeds EA ier tas Feith Tey ett ent ' ran rte oan 
Bs H 
EES sae oy = iy? PRBS See com eR a Te ie Lee Ber : Hae SELEY Hh es A : tm es . » 
NORM te 4 ee aes ert + ae Ph cae StU «Un pa CA terete ses aye \igeee ny SOULS siete TR Or Cec I ’ . ey Cate eee, ee 
Pot aaa ake +48 Lee Ae i} mie ee Mita tart . eae 5 Dt is I aN ar) a, er) ed ee r L : 
a Si Os Pay . ba ald ate) BY kas "i tar reat i he ale 54 Sd CL ea ¥ittrwenane | . > a 
* , Aa Seoge eth 4 A aie dn Bris PIS a etne bi, Soh ten SCD Dah " NM 68 Ps Fs P '‘ a a 
WS aust mes SEB A He rays boyy ft Ha te Pe atsn eaten * em bbege she i anee F Fi 
rn dSyasie gy em abel” be ¥, Eft WA ape tg oe Thee ae al har hd eae Ty a t= > fis e b 
BY, PEPE Da tae yobs tap. Pweeve Maria ts WS atep waa rae Ea oe ask rH hf Ae tite Me. ioe . ah OUAC ST 
eral eta td ee eer oteT eh orn ee eon se Or ear T Sra SOLE ICE Har epee P : ea : ety 
ee) ors eee * Pe hit eet ee eee A +, bps eas Palys ye wbey) 5 Puede be + ee ree ‘ © ee okey bey wal 0 Lares EOE . (Mesee ar) y*Pe 8 re) 
Kee eres Cae Soars bie A ale) rs baad DUT Ie Cece TEUDY. ty ee or | rT: > yr Tn Iie t ae Brice - , CU a 
beta ee Datla Ut att a rea ee A Phebe eta i OL Meas Ga tace Uh a teria ir tea Rar anal ater it, “Neb “iey eo Ep cages A 5 = 5 
Seton cao: vs Sree AAT hea daee SU reek A Pear at Dea ee Cert cd tame ne Raa ee Maron ei igi Be 7 ‘ + 
eb IS we Lead i ay nd aes bd b.) % : bed 3 A A . i bh Ld LJ ‘9 
wk. aR eae ed 1? RR eate Raae 18 at “3 een 5 Ee ey ce Hae eS ON me pb 7 4h Ha pe Oy. WEF SEG Da H 
CF e . J * Ot Fi 3 on 
es a oan eH “4 A ba : Led Se tise) AEH Bias 2g Seed Petar ria etd Be en Ty a Tee a a Pa Se ak PTY Sarva teh y 
2 ae nee ei tocerg bah ag atect pants? cfd Ree eS a aber Te at) yaetterge te eaees 
f a sf gt "YF oe My ye ee ; os “a bt ra ae 
"¥ ater > MPytiee Te Le Ae Pa Fo | a CCU rye U 
aerate ast beat Sa ae Os tr aan 1 hel ee a He A” Pande “P 
Pet Laon as ta Wy M4 bas: F a wy rere yt a ory ae rie oer an 
“rsa, ONE es Fe PY veY Ait: 2 ite ae ‘ i ROA Bey 
5 a hag nga es Dee re 14 ¢ i te Phos J 
da nee FT] RRM Ps i oh eh Ras Pt SY 
Forni FHA y . 


bes tree 






























































Le Lot ee ee | 
He cays hve & 
fl | Ce ee CP a ee Bt 
ear I ie Se tie Ca 
! HRCA unr ant her ae 
Ge Sake anuepe ince ete bs eh Lt ae 13 pier rer SLs sus cenrs Ba Genel ‘ 
ee om iM dats CCRT s i a x NOY Rhee Pith Oat oe Aare obey Eh gre Stet tee be 
p C8 id eee WR hy aT) Ta fetus tly Lee 24 Laman Sears grit aty Wa MEN pele! pa ie oe 
A ~ 2 Meee) Ct Hen perio Fadnga se eh WAP Sse - tf 
a ey his ee Uk a4 BY aes Aye Mat dee WE That CEP Ne 
ide dips Se ' ) Bl eed 
ee Ets Ae iy 
She ol al Sieahaptie) Tot et 9 r 






rf +4 Doyryte ween, Dee 
rk ee Rao Pal 


SOU ey tate og 













* ose 






™ 4 vonsuracy Shy 


Ae shay Cae ere ee pale Deter or 


ea haa arr aes yt 
STURN Ey ree Ye fey rr? 





$e ePN S Per Bere ara Vege DF per 
edie at | Oe I ee er er 
i 





Al 
Sout aeatnle 


ae Te a> is SoS PaPinctate pod 











































































































































































foes 
a Pen ad dagony. PL Card 
Hea a Oe | OPT Bye bo besa ye 
, eee oer a he ries 
Siete BY 2 he ors fe Serene eet ty Rear machetes ; aA aR Hea ah no i ae ae , er - 
we Dae ee: SOT oF | Pe aD HER a i , PEO ee ar ree nn Lat A SI? redid hh Sirs ia 
RDA ae a Ee Sa Rt De eer ee Reece SOU Ua nes on Mure ond aed rhea ee Tee | 
i Fy rede ATE PL 2 ea SIR ET eE wrt: net $ 5 aa oe 5 F) ray gee a aes 4 
ea ee aS lala CO Beh Sr ange Soa Th Gs Ae AES arty cro eat = cater all te ho SEP mha%ene oom oe ee Syste tell lg “4 
3 5% be ih ae eh So at \ t Perse ty Peoria pare tnn 4% ‘ a Aa om on pace A eae 
ist) ea Sou ae FY ede Herd Ae ee Sarg hce HiT ad “a4 MPA ate ‘ d r ‘er 
a Kets. ne alt $ aa Per Ts Hs Pret Fy 
an Shon aya UA Ot Ses Be Rts ben Me be Sed nk Pd A Piet A Ct Le villains praise : oat orp 
ey ‘ rd Uh Bd ae ea eae Me aN Ste a Bp tett ves hed ie Ae Hep 4 GET owe Teoany 4 An 
BEER renee atti Ad Ache ee ar ati ae cere ee te ACTA ns atts et! 
ts) ae at i STR eo at i a Ue rece rt i Bs ae . f a a HeLa 
SEE Le ee iS eit, ay phate ues cat pea Sa ie ite a cares SW ah ea WOH oe tae t et ang ently oe Ar n aro 
axe Ka spe Te eee tan pale aie WATER Mare Hoe ae ae DUA RECUR PL RS Poi tere eo enene 
a an oP Ree mitt fay) hte rs pay site ih ry Q1 cee a . ry . HW op atarde A . aire a AL - | Le bl 
be siabh te Matar raed er net Ai Hie at Sy a Pe Py Sa eer ee rg oe are wrt hepad eet LN ada IY) Seer YT aL tela ra Pe we JA ats Che Lec ‘ es 
Fonwesy an aeatued apa ts a CAPE Ryn sete pr: eae hits chetans Thera a teea i a ea re ae ae Sl ACCT FPR Poe aan 4 F ee ante ati dareve a 
oaraek Gis ght efit aos AURA Hdd et Sled Aa Te te a i ree Cie Ais be NaS tery} DU PEeraT TRY erry, area SS a eee Le 
Se tate ah 0s ae Pett et tthe i, Air Spe Satie od Th pia kr eta bie oh sree " pad) Ald Sind SH alt ie LATOR ry rd eh? ee er ee 
tiem Beet e nary tl ys, ty Coy a $44 bf Blafete pre Mi We - Tait Wycee et 
ple ee ees Aoi te ad ae en : aL Tree NL Tate he Res ea “rt Men de sles ae 
Lian pte lan Is are ye raked Lies rogh DUP grr he by tet Lh te i 
aE eee sees Cate g te i 28 raed rive auhesispten Se 
i ( Sh laid oa tol TTT TY be On 
4 A 
a Pale 2 tetas Ted it 


TGs rea oe ra 
, Bes 


















rt “BOe at yey ie sep faite Th 
Le ut haere 

Oe Nias x Tog 

FA WP Pi ee Y 









de 
























































































fe "ehtars rr 4 
rman verge Be 5 Fe} La 
33s er At NO e' "974. iybs ang % er eras oH 
Py phan 45 imag y oad Sate CO Are th: Ce oe er ee Ln ae! 
2 = re haa sy tee pees rah) see re Le mre i ech Meera rr eae ET o 
Aer > ae oh t ie fy M4 2 5 Pere ees ec ieee 
aon ee sgn ahr, ru +t Patients Ch Ei wnt Mea a Tot en ; mets Ciara yay ae oe chor La ated oet re a a tae 
v Liking tel ° cy rh i De had , Lite la] is cra ery er wy Py . Tan 4 
‘1 i Deberete tT Mien a :! are renga s Gi eee yy arakes bi ta 1 yeast PF a mye aris ue Pa 5 ods Pa 
eh bh ta aad eee te oe AES te a art CoS Pe i} at a Ps tet 
nal tt} Mbt Wa Fs Pere ye Haat }4 SAS On ie a 
4 rs " Aes ea EN ‘ 
Hee peta ks nie Ware te 











eee 
. x a 24h WT Oe 


ts ci) 


















Ay bee si 
Para a oP td ee 






































































yy 
; ae 
2 7 Par Pan Yi our A bred ype nee ks 
vy ! ean er by ae aa a etl f parte ae a Sar hdde HLL one Aas it dey aber api Pe ae 
2 rr it Od as ¥ Ba if ney te ieee eg a) a eee ee oan at Pee eke 8, » 
ai sme mee rtd 3d tas eS Ya gt TUR or ford Mure Bere 
fot Se St bd> a tar AI AL eo fey H 
. Pt bold Lat) ‘ Area a bar, n 
Fae as Pad by tt 0 " 


ape Pert yey 


1 Le 
iia ashe ay He Sg San 


a ark ok ror eres weer yd oer ‘ 


© 49t Fs 
a 









eee: res 
ete fs Showy 





ne ee 
me ae Ca ry 





herd tis 
re wees rn 
7 











ees, 





4 auth axe 
aS tn i} rer mS 

Bt APO eAESe Nit yretedy Ha we re ey Nalary 
ee San Loa Oh Oh es Merssgrtauthas 
poe ee eee cy reas arte eae? rat a 
a oe Rihana oe iT yids y oes Oar wit aee¥, th bat ies Eres 
Prony "ss hd olla att ett aan rr) es sh * ss 

* ee tLe Ree VSS my Sate ar et ove 
































ah & 
PVE PAG huey ae ne 
an A) a eet 
bt Pons) WU poe Crd 

od ed | oe Ons 

ta I 





















te td tt 
RATAN Rete N TTS 


ey at 
eat TS dott ute oy ease} rer 
















pad bags aye % 



































































































































































‘ * ay r 
vee neta ny rant i: 1 @ae ase 
a ae eget ete eo : , tae aa Cea Pepe ae ee 
ace ne es . oh 
ME os pasate at sete Pe U he we aes $ Weedk ya, an 53 Hanes O aL ee Bi war 5 : r Pat a 
Peper bet seat Ogi ats ay Se Patan a Lt wee Gs6F nL Sin fp SC CAE Ur err oar) 
aie tte ido A Perc CNS Reet ae Py A 0 Be Sets bibaat +1 ee *) SHG ue al be Ary t any en RT err re x a ei ge 
¥ oe ei b Mabey F A A 6 
at ree oad ch Stee Pear et ena ee TE e ae on See ee a ee os Bega k 
tw tA Sed Fj i ae rer mi Ce he a oe *ute ry ry 
ae nea lds Oats 4 Fpl Rees bi he ' + eo crehes Mate ot ie ’ ir lag Leary *, ‘ 
os Pi ident os bea ty om PO Per ere ere we ig etl ots ae mye ame Ug ren A Ny O ry *¢ el 
4 he ks reser poe iY at fe hy ees ele Eee) 4 ae Tmt - a! é 
rat eae te Fn Peete eS AM Ae g? Terese A fare a 5 ar eh they a «ork Wh ore Ue CAL ae 4 ee CL Ga 2 ebws 
peat aya te re Bits ws aa bh LS a Ae neers ine Lares Het 4 Vode f Lien Le ers B rs Lyd are y Pay ‘I 
a esr eran ake tt ae Sah era hs ae Srl n eas es 1 i Hn id] ¢ q™ Ha) i | rT a 74 oa Toys 
is EN RS fens ia PLS AR DRA RE ae MTSE SE AR 5 Ree te Bee 
‘ Sint at iene oy mera Sart oy A 
Pees ie He aa GORE OTT EPS TT Hee Ais Pee het RETR et, 
“1 ORR A db HOWE Se Pn ists May Fee! ~ ia PAT 
i 7s wiazss yt Sitar Der Ge LP Nua Pie ee TS ( 
YEG Dik ry EN 
Lins b abla | Pereltebs 
Pere ty ee 
eed Py oh 














































ob, 
rt Tata pe ete Wie 
A * .. i ra nile 5 ee > ; G a a Pd ee 
+ AD 2 bd Vy 
id Raitee et tat als Ha ; au MP us etSkes 

paerent a He tata vie aH Rat elt roar 

Ud — 
Stites 








Lo] yes bo 
a RMS | ih : pee rs 
Banat eaten Paella 

esti 


tis ve 
Pay eae 
ST ‘a 
‘ whee 4 

























































sl at ia i 
t b 1 oa eee # ars Uy 4 4 
ve PAE Menara Patch pepe ‘yin Se vie ie he 4 ne 
H Ad blared be gt f Mahi he ie ia 
Bears ed 4 


- dee 
a 









ARETE ot i 
4 fyen- fy a 
rs ay etangreceaee serie k 
See BLS Oe sae ee 








enn Fas TREN at 












ret 
i) 
oe yer Lal a Piers Pe Rt e y ah 
a ee ND Ae ed 4 aT 





Sneha wy 
BEd 
tty 





Pia Cr 
hn RON 












ae 



















hed a) ch ol hon 

Cte at, aA 
eerrg Aeon H 
Lad to 


Neos unis 
Wat B ihe eto 
ay ky Pah 



















it 
WA Ca 
Y 







ba pr oe Ww) 

Fae Te ER 
perere! ar hie 

eres y: Cy Rk ES 


F Noeram 





ea 





Pes 























, sr 

| 

b : tT on a 

h, yt oe ae -s ie ele 
yy 0 

- Fert « tv 

+ 


